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Properties of linear passive media are investigated by using a phenome- 
nological approach. Properties of fields that can propagate in a passive 
medium are postulated and from this properties 6f the constitutive relation- 
ship are deduced. A necessary positive real condition on the constitutive 
relationship is found and some of its implications are considered. Also, the 
causality condition which is necessary for realizable media is considered. 

Next a general formulation of the spectrum of characteristic waves in 
lossless linear passive media is made. Because of an orthogonality condition 
for the characteristic waves of the medium, the fields due to an arbitrary 
source can be separated into components parallel to the characteristic waves. 

The components of the source field are dependent only upon the portion of the 
source parallel to their characteristic field and to their own sheet (s) of 
the dispersion surface. The theory is then applied to two particular problems, 
electric dipoles in a general time-dispersive uniaxial medium and in an iso- 
tropic compressible plasma. Finally, the radiation field of an arbitrary source 
in a lossless linear passive medium is investigated using the spectral decom- 
position of the fields. By normalizing the length of the Total Poynting vector 
(electromagnetic plus medium) to unity for each characteristic field, a concise 
and physically interpretable expression for the source fields is obtained. 

These results are then applied to an anisotropic compressible plasma and to a 


magneto-ionic plasma. 
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1. INTRODUCTION 


One of the most simplifying and useful properties of Maxwell* s equations 
is the linearity property. However, the linearity property in the past has not 
been exploited to its fullest. Of course Maxwell’s equations will not be linear 
if the medium is not linear. Thus, in all that follows a linear constitutive 
relationship will be assumed. 

To the author’s knowledge, the most general linear constitutive relation- 
ship that Maxwell’s equations will allow has never been used. There is, however, 
considerable motivation to formulate such a theory. In most cases the permit- 
tivity, e, and the permeability, p., are considered to be scalar constants. This 
was probably the original formulation of the constitutive relationship. Then 
these scalars were found to be functions of the applied frequency which gave rise 
to time-dispersive media. From this, it is found that the group velocity and 
phase velocity are not necessarily equal but are in the same direction. In crys- 
tal optics the permittivity constitutive relationship is given by a diagonal 
matrix, e, while the permeability, pi, is still a scalar. Two such media that 
occur in nature are uniaxial and biaxial crystals. More recently a considerable 
amount of the literature has been devoted to magnetoplasmas which occur in the 
ionosphere and can be created artificially in the laboratory. Magnetoplasmas 
are characterized by a matrix permittivity and a scalar permeability. They may 
also be space and time dispersive, i.e., the permittivity matrix a function of 
the Fourier wave vector, l£, and the Fourier time number, u>, respectively. Also, 
there exist materials such as ferrites that are characterized by a scalar per- 
mittivity and a matrix permeability. Further, uniformly moving media has an 
even more complicated constitutive relationship. In a moving medium, the 
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electric flux density, D, is linearly related to both the electric intensity, 

E, and the magnetic intensity, H. A similar dependence results for the magne- 
tic flux density, B. These are a few examples of well-known media of nature; 
however, the possibility of creating artificial dielectrics gives impetus to 
investigating the general properties of linear passive media. Moreover, the 
synthesis of media with given space and time-dispersive characteristics would 
be highly desirable. 

A large number of workers in the area of magnetoplasmas attempt to derive 
a macroscopic constitutive relationship of the medium by formulating a micro- 
scopic model that is describable by dynamical equations and/or probability dis- 
tributions. Thus, they attempt to deduce the macroscopic properties of the 
medium, and hence the fields that can propagate in them from a postulated 
microscopic model. This is not the only method of deducing properties of the 
medium. Another very effective method is to postulate properties of fields 
that can propagate in the medium and from this deduce the properties of the 
constitutive relationship of the medium. This phenomonological approach will 
be the method used to investigate the properties of passive media. 

Since Maxwell’s equations are linear, the spectral theory and represen- 
tation of the operators that compose Maxwell’s equations for a lossless medium 
will be useful c With this the spectral representation of both the source-free 
and the source fields can be derived. In this way, the specific role played 
by the sheets of the dispersion surface and their associated characteristic 
fields is easily seen. Also, the manner in which the source excites the spec- 
trum of characteristic fields to form the composite source field can be deter- 


mined. 
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In Section 4, by using the spectral representation, the exact expression 
for the fields of a dipole in a general time-dispersive uniaxial medium will be 
derived. These expressions are then compared with Clemmow’s^ fields in an 
uniaxial medium, which he obtained by a scaling procedure. Section 6 is con- 
cerned with the radiation field of an arbitrary antenna in a lossless medium. 

The purpose of this thesis is to exploit the linearity of Maxwell 1 s 
equations and thereby derive general relationships for linear passive media 
and for the fields that can propagate in such media. Further, the usefulness 
of the general relationships are demonstrated by specific examples. 


2. PROPERTIES OF LINEAR PASSIVE MEDIA 


2.1 Maxwell's Equation s as_ S i x- V e c t or s 

Most commonly. Maxwell’s equations in point form are expressed as two 
separate vector equations, Ampere’s Law and Faraday s Law. The vectors in- 
volved have three components. However, this separation is by no means unique 
nor necessary. With equal ease, one can express Maxwell’s equations as six- 
vector equations involving vectors with single components only. 

Likewise Maxwell’s equations may be expressed as a single vector equation 
involving vectors with six components. This is the formulation that will be 
introduced presently and will prove to be of an ideal form for many purposes, 
particularly for conservation of electromagnetic energy. Define the partial 
differential operator O and the six-vectors } 3^ and ^ as follows: 


& 


0 Vx 
-Vx 0 


( 2 . 1 ) 


3 



( 2 . 2 ) 


3 




(2.3) 


S’ 



(2.4) 


The subscript "f ' is meant to signify flux since the units of involves 

per square meter. With these definitions, Maxwell's equations are simply 
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2 # 2 Definition of Constitutive Relationship 

As easily seen, this vector equation has more unknowns than its order 
and hence is indeterminate. The constitutive relationship yields the remaining 
equations of the system. Essentially, the consitutive relationship is a 
statement of the influence of the material medium on a wave propagating through 
it as compared to the wave propagating through a vacuum. As such, the source 
is not involved in the constitutive relationship. One would like to formulate 
a mathematical expression for the constitutive relationship that would encompass 
the general class of linear media. In the real space-time domain, say R, such 
a constitutive relationship must be a matrix convolution operator, denoted ?/° . 
Then the general linear constitutive relationship is 

vr * 7 < 2 . 6 > 

where the convolution is with respect to all four variables of space-time. It 
can be verified that this expression is capable of accounting for the scalar, 
uniaxial, biaxial, magneto-ionic, ferrite, and uniformly moving media that were 
previously mentioned. Furthermore, the general linear constitutive relationship 
applies to transient fields as well as steady-state fields whether sinusoidal 
or otherwise. 

In light of the complexity of the convolution linear operator, which in- 
volves integrals of the "operand," $ , at times it is preferrable to deal with 
the system of equations in Fourier transform space, tT . Throughout the fol- 
lowing, the definition of the Fourier transform and its inverse transform will 
be, respectively, 

F(lT.w)= ?<M)d 3 rd. 

-CO 


(2.7) 



6 


and 


oo 


y(T^)={2if) A JJJ V 

-OQ 


j(wt k > cu)d 3 kdOJ 


( 2 . 8 ) 


where the transform variables are oo and k. Let the Fourier transform of 
^ y , and 0 be, respectively, F , C, U, and 0. Then in space the 
system of Maxwell’s equations and the constitutive relationship are 


OF = OJF f + C 


(2.9) 


and 


F f = UF 


( 2 . 10 ) 


What is gained by such a transformation is the following. The operators are 

simply square matrices of order six with elements which are functions of the 

transform variables co and k. Also make the distinction between space-disper- 

2 

sive and time-dispersive media as Allis does for magneto-ionic media; the 
constitutive matrix U as a function of k and co implies that the medium is space- 
dispersive and time-dispersive, respectively. In truth the names space and 
time-dispersive media are misnomers since dispersion is a property of the whole 
system and its dependence cannot be separated among the space and time varia- 
bles. A non-dispersive medium is one whose dispersion surface is hypercones 
with their vertices passing through the origin in space 4r. Thus, it is pos- 
sible for a space and/or time-dispersive media, as defined above, to be non- 
dispersive. 

2.3 Waves in Passive Media 

The general class of linear media is larger than necessary for most pur- 
poses. In fact the homogeneous media that commonly occur in distributed 
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systems is passive. (A passive medium is defined to be a medium in which both 
the electromagnetic energy density stored and the power density dissipated for 
all fields are non-negat ive. ) For this reason, in all that follows we will 
restrict our attention to passive media. First, however, it is necessary to 
recognize the properties of passive media in order that the two classes may be 
separated. 

For either passive or active media, the totality of fields (waves) that 
can propagate through it uniquely characterizes the media. Even certain sub- 
sets (complete sets-subsets of the propagating waves that span the totality of 
fields for the medium) of the totality of fields may suffice to characterize 
the media. An example of a complete subset is the set of plane waves. Even 
so, all plane waves are not admissable in passive media. Thus, we would like 
to determine the properties of the waves that can propagate through passive 
media and from this determine certain characteristics of the constitutive rela- 
tionship. Neither Maxwell's equations nor the constitutive relationship will 
aid in determining the class of fields we seek. The remaining equation that is 
essential to determining whether or not a wave is propagating in a passive or 
active media is the conservation of energy. Multiplying Equation (2.5) on the 
left by 7 and transposing terms gives the equation of conservation of electro- 
magnetic energy. 


or 


T T&7-- (2.1D 


dvr 

d-C 



o> r 


( 2 . 12 ) 


where the superscript T signifies the transpose. Identification of the cor 
responding terms results in 
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; rate of change of electromagnetic energy density 
supplied by the sources. 


-7 T &7--V'0> e 



dw- & 

<9i + / ? de 


; divergence of the electromagnetic Poynting vector. 


; rate of change of stored electromagnetic energy 
density plus the dissipated power density. 


A further identification of the rate of change of stored electromagnetic energy 
density and the dissipated power density can be made provided the operator \f 
is split into the sum of two operators. In the transform space, , it can be 
shown that the rate of change of stored energy is associated with the Hermit ian 
part of the operator U, while the dissipated power density is associated with 
the skew -Hermit ian part. Thus, it is only necessary to transform the Hermitian 
and skew-Hermit ian operators back to the space-time domain, R, in order to ob- 
tain the separation there, i.e.. 


and 



j/v,t)*r T (-7 



( 2 . 13 ) 





( 2 . 14 ) 


Then we have 


and 


7\V % * f) - 
7 r (¥ a * 7) = 


( 2 . 15 ) 


( 2 . 16 ) 


Now the separation into passive or active media can be made by determining 
whether or not both the stored and dissipated electromagnetic energy densities 
of all possible fields in the media are non-negative. This will be an essen- 
tial factor in proving a necessary condition "positive real" theorem for the 
constitutive relationship of passive media that will be considered later. 

At this point let us determine a property that all plane waves in a pas- 
sive medium must obey. For a source-free passive medium, the net energy flow 
into a closed region must be non-negative. Since the closed region is arbi- 
trary, this implies that the integral with respect to time of minus the diver- 
gence of the Poynting vector for any field is non-negative, i.e., 

-/V^dt>0. 


Also, since the field is arbitrary, the condition should apply to a plane wave 

cl the type 

3 - 2 Re (2.17) 


where 


3J = 3£e ,So, - 7 °'° 


( 2 . 18 ) 


?o is an arbitrary complex vector field which is independent of space and time, 

s o s< V*j«o, 

-v>. ■ 


( 2 . 19 ) 
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where 


Therefore , 


or 


U ° l-Yo* 0j 

-VV e = Re[^0 0 ^+^ T O o ^j 


-V ^ = 2 Re 


( 2 . 21 ) 


( 2 . 22 ) 


'f o 2(^ c t-a 0 7) T 2(s 0 t-y c -T)l 
^0°0^0 e + ^ O 0 0 ^c e J (2.23) 


Let <j > 0 and time, t , be much larger than any characteristic time of the 
system; then the above inequality is approximately given by, 

f^Sj T /p Mi 2(0i)t- ct 0 l') [" I ✓v I ^ T n 't 2j(w 0 t-k 0 

1 V(r e dt~e Re ^ %0, J 0 + -^of 0 0 c GF o e ' J >0 (2.24) 

2(Cr o t -d-c^r) ^ 

Since e >0 


then 


Re 




1 </- < 


2iK. t " k a ,rJ 


■ r> I 

> 


0 


(2.25) 


The second term of Equation (2.25) may take both positive and negative values; 

s 

therefore, the first term must be non-negat ive r In other words, 


But Ro 


[^O^o] 


Ke 


-4Qf 0 P e where 


[•?to 0 3,] > o 

= yRerE 0 XH*J 


(2.26) 


is the real part of the 
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complex Poynting vector. Now we have the desired result that for Re s >0 

o — 



then a Q > 0. The result is as intuitively as one might expect; in a pas- 
sive source-free medium the fields decrease in the direction of the average 
power flow. 

In addition, since the conservation of total energy is also applicable, 
one has 


dur r 

” 37 “ 


+ ^ 1< ^ > T 


a^r 

at 


(2.27) 


Therefore, the difference between conservation of total energy and the conser- 
vation of electromagnetic energy is also a conservation law. That is, 


dvr^ 

at 


+ 


^dm + V T <P m = 


at 


(2.28) 


where 


- ur e 


(2.29) 


P dm ^dT fide 


(2.30) 


(P m = (Pj ~ (? e (2.31) 

and 

im = i T - i-e (2.32) 


The subscript m is to denote medium energy density, medium power flux, and so 
forth. 

Naturally all of the conservation laws may be expressed in transform space, 

T ; however, it is sufficient to only deal with the transform of the conservation 
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of electromagnetic energy. Then 


g. -p 

~ ~ v I /T) 

sits ' q + ^ cie ~ ' (s Q - s ^ ( 


(2.33) 


where the tilde indicates the transform of the quantity. 

2,4 Lossless Property of Media 

Consider a harmonic electromagnetic field of the form 


- Re 


j ( ol> t k ° r ) 

v r\ S 


(2.34) 


with 


7 f = Re 


j uj U ( k , co ) y Q e 


j ( cu t - k ; r ) 


(2.35) 


By averaging the quantity J ^ with respect to time and equating this steady 
change in energy to zero, we can arrive at the condition for lossless media. 

The terms involving e — will have a zero average, leaving 


y\ ( jwy(it,cu))y 0 + yl ( j^y(if,cu))V c 


= o 


(2. 36) 


or 


joj ^ 


yCiT .uj)- y + Ck\aj) 


^ 0 r o 


(2.37) 


But is arbitrary and hence the constitutive matrix is Hermit ian for real 

go and k , i . e . , 


U (]T,aj) = y + ( k , oj) 


(2. 38) 


2.5 "Energy" Condition and Group Velocity 

Consider a lossless medium. Hence, for real co and k we have, 
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U(k ,cl>) = U + (k ,u)) 

and 

0(k) = -0 + (k) 

Therefore, the matrix operator 7^(k,co) defined as, 

7T){k ,(jj) = 0 (k ) -jCL»y(k,cj) 

is skew-Hermitian. A source-free solution to Maxwell’s equations, F, 
r exists when 

77jCK ,oi) F CK.cj) =0 

or 

(^F) + =F + Z7 + = -F + ^ = 0 

Allow small perturbations to occur in oa,k and the medium such that 

Cut, = Cl) +8 CJ 

and 

k, = k ♦ 8k 

Then a new field solution F^k^c^) may be obtained for the perturbed 
such that, 


( 2 . 39 ) 

( 2 . 40 ) 

( 2 . 41 ) 
in space 

( 2 . 42 ) 

( 2 . 43 ) 

( 2 . 44 ) 

( 2 . 45 ) 

system 


77|, (k, ,U), )F,(k,,a»,) = 0 


( 2 . 46 ) 
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7T| | is not required to be skew-Hermit ian. Because the perturbations are small, 
however, we can write 


7^ ,(1^ , 0J| ) = 77JHT, OJ ) H-Scu ^(|T, C4 J) + W(T,cli) + 8^(Tw) 


(2.47) 


where 8 Jfji U OJ ) represents any perturbation of the medium not due to a change 
— 

in co or k. Therefore, by using Equations (2.46) and (2.47) one obtains. 




F, = 0 


(2.48) 


However, from Equation (2.43) it is seen that 


-t, 


F { rf]F { --0 


(2.49) 


and hence 


’[ Sw ^r^ +sT ' s 'k n l +5 ' r '| 


f=o 


(2.50) 


Since all of the terms in the brackets are small, the difference between and 
F may be neglected, yielding, 




F : 0 


(2.51) 


Now let us consider each term of Equation (2,51). 

F+ ^^ F M F ^(wU)F = -j4W (2.52) 


W is a generalization of what is commonly called the average electromagnetic 
energy density due to a harmonic field. This is discussed further in Section 


15 


2.9. From the second term we have 

F + [v k ?rj]F=4j[p e *p m ] 

where , 

P e * - -y iF f [V k o]F = l/2Re[EXH*j 

= electromagnetic Poynting vector for a harmonic field 

and , 

W p [^y] ^ 

= medium Poynting vector for a harmonic field. 

The third term reduces to 

F^jsT^ j F = - j u) F + 1$ U j F 

Putting the component terms back into Equation (2.51) gives, 

+ Pm ) = § ^ W + i- uj F 

Now for an unperturbed medium and we see that the group velocity is 

given by, 

(2. 58, 

Also when 8 OJ = 0 j i # e. , when the refractive index surface is given for a 

fixed frequency the group velocity and the total power flux vector P = P +P 

Tern 

are normal to the dispersion surface. It should be emphasized that the only 


(2,53) 


(2.54) 


(2.55) 


(2.56) 


stipulation made upon the medium is that it be lossless. 
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2,6 Onsager's Property 

For the sake of completeness let us mention Onsager's principle as applied 
to the constitutive relationship. If the sixth— order matrix W is 


w .ri °i 

= -ij 


(2.59) 


the Onsager's property can be expressed as 


U T (k,cd,H 0 )= WU(k,CL» ,-H 0 )W 


(2.60) 


4 5 6 

For a description of the details refer to Onsager, Meixner, Casimir, or 
7 

De Groot . 

2,7 Real Property of Media 

Since both the field strength $ and the flux field are real valued 

vector fields, the definition of the transform and Equation (2.6) implies that 
there exists a real property of the medium, namely, 


y*(s,y) = y(s*, y*) 


(2.61) 


where the star, *, indicates the complex conjugate, and s = <j + joj and ■y ’ = j~k^ 
More explicitly, $ and real implies that 


F (s,y) = F(s , y ) 


(2.62) 


and 


Fj (s,y)= F f (s*y*) 


(2.63) 
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respectively. But the constitutive relationship implies 

F f (s,x) = u (s,x)F (s,x) (2.60 

Therefore, 

F*(S,X ) = U*(s,X)F* (s,X ) ( 2 . 65 ) 

and 

f f (s*,x*) = y (s*,x*)F*(s ,x ) ( 2 . 66 ) 

Subtracting Equation (2.66) from (2.65) and using Equations (2.62) and (2.63) 
results in 

°= [y*is,n-u(sV*)j F*( S ,n (2. 67 ) 

The vector F ($ f X ) is arbitrary, however. Therefore, 

U*(s,X) =uls* X*) (2.68) 

An immediate consequence is that if the constitutive matrix U is independent 
of both s and then it is a real matrix. This property does not seem to 
have been fully exploited in the literature. 

2.8 Positive Real Property 

It is quite well known that the necessary and sufficient conditions for a 
function to be the driving point impedance of a linear passive network is that 
the function be a positive real function as originally defined by Brune 8 in 
1931. At that time Foster had already presented a method for synthesizing 
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certain types of one port networks. This greatly reduced the difficulty of 
the sufficiency proof. The significance, however, of Brune ’ s positive real 
theorem was that it established a logical foundation for and an impetus to the 
development of network synthesis. Also, it has been established that the 
impedance and admittance matrices of passive networks are positive real 
matrices. With this as a background, it would be desirable to establish 
similar necessary and sufficient conditions on the constitutive matrix or 
perhaps some function of the constitutive matrix for passive media. The de- 
gree of difficulty of such a proof over that of Brune ’ s for networks will be 
enhanced for at least two reasons : (1) Very little is known about synthesis 

of media, particularly the synthesis of dispersive properties; (2) The proof 
will include the wave vector k as a parameter and thus will require extra con- 
sideration and conditions upon it. Heur ist ically , one might expect that in 
the limit as the wave vector k approaches zero or as the wavelength approaches 
infinity a correspondence exists between the properties for passive media and 
passive networks. 

Let us now establish the necessary "positive real 9 ’ property of passive 
media as a theorem. 

Theorem 2.1: The function of the constitutive matrix Z(s,Y) = sU(s,^) is a 
"positive real" matrix for a linear passive medium^, where s = Cr + jco and Y = jk 
are the complex transform time number and space vector, respectively; i.e., 

(1) Z*(s,y) = Z(s* y*) 

( 2 ) Re s >0 implies Re jP^Z^ > 0 , 

Note that the restriction of Y = jk where k is real is necessary for the convo- 
lution operation in space-time which gives to exist. 


4 
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Proof: Part (1) follows immediately from Section 2.7 and the definition 

Z(s,y). That is, 

* * **“ -^ ♦ * ♦ * .«»<* 

z (s,y)=s u (s,n=s u (s ,y )=z(s ,r ) 


Recall that 


dK 

dt 



(2.69) 


Both2^(f) and ^t7,j e lt) arpcpositive definite or semidefinite for fields in 
passive media; therefore, 

t t 

2 ^> de dt = y*^ T ^fdt + 2(T e (t 0 ) >0 (2.70) 

*o <o 

This is true for any field; thus, it should apply to any particular field. Let 
a particular field be ? = 2Rc 7| where J^j = ® ® and ^ is a 

complex field vector independent of space and time. The flux field is then 


or 




(2.71) 


Ff = UF 


(2.72) 


Now 


F = 7 0 $(s-s 0 )S(x-r 0 )+ ;f*&(s-s*) 8 (y-y*) (2.73) 


and 


F f = u(s,xi 7 0 8 (s-s 0 )S(x-y 0 )* u(s,y) 5 f ; 8 (s-so) 8 (y-y;)( 2 . 74 ) 


Therefore, 


^f = y(s 0 ,y o } A> e + u(s 01 x 0 ) 



* 

(s 0 t-y 0 

e 


r) 

(2.75) 
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However, because of the real property, = U*(So»)^) > ^ ux 


field and its time derivative then becomes, respectively, 


^ f = 2Re[u(s 0 ,^)jf 0 e lso, ' ) '° r| 


(2.76) 


4 =2Re[2(s 0 ,T 0 )3! 0 e ISO '' 70 ' 7 ’' 


Now it is easily shown that 


^ = 2Re[^2(s 0 ,y) 3,+^Z(s 0 , y 0 ) 3- j 


(2.77) 


(2.78) 


= 2 Re [? 0 Z (s 0 ,!j) 7 0 e^'+jf z^^e 21 ^' 


(2.79) 


Using the expression given in Equation (2.79) as the integrand of Equation 
(2.70) and integrating results in, 


R e^^o z(s o'C 1 -^o e2<T ° l+ TT 4> Z(s °’ ) ° , ^° e 0 ° ] 

- R e[^7 0 T Zls 0 ,^)^e 2<r « , ° + l ! r X T ZU 0 ,^)? oe 2<5 o l o- y o r| ] 


(2.80) 


+ 2^(t 0 )>0 

For a >0 and for t large enough, the values at t i.e., exp (2 o t ) and 
o ? o o o 

exp (2 s t ) , are small compared to those at t. Also the electromagnetic 
o o 7 

energy density at t , #£< to j, is small. Then 
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and 


e2V r, ]>0 (a.8i, 


Re [oJ- ^Z(^i^)3( + JjjfJziSo.^J^e 2 ' 11 ’ 0 ' k o rl ] >0 ' (2.82) 


Since the second terir. of Equation (2.82) takes both positive and negative 
values, the first term must be non-negative, i.e. , 


Re s 0 >0 implies Re >0, V3« 0 ( 2 . 83 ) 

Therefore, as was to be proven, the matrix Z(s3 = s JU (s , Y) is a "positive 
real” matrix. 

In practice it is difficult to prove or disprove that a medium satisfies 
the positive real condition from the definition alone. For this reason K-t 
us prove a lemma that will facilitate this work. 

L emma 2.1 : Let Q(A) = X + A X and QCAp) = X + Ap X where X is an arbitrary 

vector, and A^ is the diagonalized matrix of A Then Re s > 0 implies 
Re Q(A) > 0, VX if and only if Re s > 0 implies Re Q(A D ) > 0, VX. 

Proof: Let T be a matrix whose columns are eigenvectors of A. Also, let T 
be normalized such that T T = I. Then when the vector Y of the transforma- 
tion X = TY takes all possible values, X takes all possible values and vice 
versa. But 


X i AX=(TY) + A(TY) = Y + (T + AT) y = y + a d y 


VX,Y 


(2.84) 
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Therefore, 


Re s >0=>ReQ(A)>0, VX if and only if Re s >0 =*>R e Q(AJ > 0 , VY (2,85) 


Thus, to prove the positive part of the ^positive real” 1 ' theorem, it is 
sufficient to show that for Re s > 0 the real parts of the eigenvalues of 

Z(s,k) are greater than or equal to zero c 

o 

2.9 Implications of the Positive Real Property 

Several implications are a direct result of the positive real property. 

For instance, it is possible to derive a result analogous to Foster's reactance 
theorem. Before doing so, however, let us consider the average electric energy 
density in a lossless medium, It is quit© commonly stated that the average 
electric energy density of a harmonic field at frequency w is 


<vaE + -$£ i E 


( 2 . 86 ) 


This equation is only an approximation, which to be sure is a better approxi- 
mation than (1/2)E^*€ E But the energy density is dependent upon the entire 
history of the fields in the medium. Therefore, one would expect that the 
energy density at frequency m be a function of all of the derivatives with 
respect to u, depending on the manner the amplitude rises from zero. It is not 
even intuitively obvious why the approximation given by Equation (2.86) should 
be a positive number. The following theorem will prove it to be the case for 
lossless passive medium 


Theorem 2.2; For a lossless passive medium, 


F*j£; F > 0 V F 


(2.87) 
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Proof: Let Q(s,k) be defined by Q(s,k) = F Z(s,k) F. For s = jw, Q(s,k) 

is imaginary since U(co,k) is Hermitian. Expand Q(Sjk) in a Taylor series 

about a point s = joa and evaluate it at the point in the right half plane 

s = s , 
o 




. (so-i"> 

s=ju> 


( 2 . 88 ) 


As s approaches jgj, the first term of the series will be the predominant 
term. Define the following terms, 


a = org Q(s 0 ,lT)-Q( 


(2.89) 


P ' o f 9 



(fe.90) 


X = arg 


9Q(s,lT) 

3 s 


s = jw 




(2.91) 


Therefore, in the limit Equation (2,88) requires. 


AmCt=Ai™P+Y (2.92) 

The "positive real" condition requires that jaj < 7T/2 for ||3| < 7T/2. 

Therefore, for Equation (2.92) to be satisfied, Y must be zero. Hence, 


<9Q(S, k ) 

d 1 



(2.93) 


or 


JL- F + ((jU((J,'k))F] > 0 (2.94) 
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But since F is arbitrary, we have the final result 

F+ [a?r( (j y (w ' k) )] F - 0 . vf (2.95) 

Now it is a trivial matter to see that Equation (2.86) is a special case of 
Theorem 2.2. 

A few important residue conditions that are a result of the ’’positive 
real” condition are as follows: 

Theorem 2.3: (1) There are no poles of f'*’ Z(s,~k) F in the half plane Re s > 0. 

(2) The poles of Ftz(s^lt) F on the jco axis are simple and the 
residues are positive real. 

Proof : Assume that there is a pole of order n at s = s Re s >0 The 

o* o 

Laurent expansion of Q(s) = F^*Z(s,k) F in the neighborhood of the pole s is 

o 

of the form 

co 

Q (S) =^ a i(S“S 0 )' (2.96) 

i=-n 

If s is sufficiently close to s o , the term corresponding to i = - n is predomi- 
nate. Then 

ReQ(s)^ cos|arcj(a_ p )-nqro(s-s,,)j (2.97) 

Since arg (a_ n > is independent of s and arg (s - s^) ranges from 0 to 277 for 
Re s > 0, then the dominant part of the Laurent expansion changes sign 2n times 
in the neighborhood of s ^ 9 Therefore, there exists and s such that Re s > 0 
and Re Q(s) < 0. But this contradicts the ’’positive real” theorem and thus 
there can be no poles of F^ Z(s,"k) F in the half plane Re s > 0. If s = 
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I 
I 

| then the dominant part of the Laurent expansion satisfies the positive real 

theorem when arg (a_ n > = 0 and n = 1. Therefore, the poles of F+ Z(s,k) F on 
the ju) axis are simple and the residues are positive real. 

Another theorem follows from Part (2) of Theorem 2,3, 

Theorem 2.4: The matrix of residues of Z(s,k) at any poles on the jw -axis must 

be positive definite or positive semidef inite. 

Proof: Let the vector F be real and s = ju> o be an arbitrary pole of F Z(s,k) F 
on the ja> -axis. From Part (2) of the previous theorem, such a pole must be 

-L. — ► - 

simple and the residue of F 1 Z(s,k) F be positive real. Therefore, 

Residue [f T ZU, K)f 1 = Jim ( s-jO> 0 ) F t ZU,T)F > 0 VreolF (2.98) 

L J S-t-jWo 

= F T [iirn(s-j(J 0 )Z(s7k)]F>0 VrealF 

, ^ L sH J ^ 

But lim (s-jw )Z(s,"k) is the matrix of residues of Z(s,k) at s = jw Q . There- 
s—jc^o 

fore, the theorem is complete. 

One other condition upon the constitutive relationship, which is implied 
by the "positive real" theorem, is the following: 

Theorem 2.5: The matrix of the real part of Z(s,k) must be positive definite 

or semidefinite for Re s > 0. 

Proof: Since Re [F+ Z(s,k) F] > 0 VF for Re s > 0, it will certainly be 

valid for F real. But for real F, the Re operator commutes with F, 


RejF + Z(s,T)F] = F T ReZ(s,“k) F 


(2.99) 
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Therefore y 


Re s > 0 


1 F 


implies F ReZ(s,k) 


F>0 


( 2 . 100 ) 


and the theorem is complete. 

Let us verify that the cold magneto-ionic medium satisfies the "positive 
real" conditions. For such a medium the constitutive relationship may be 
given by a matrix of the form, 


U -- 


[i oi 
lo^ij 


where 


and 


6 = 


e 1 € /7 o" 
Wo 
0 0 


( 2 . 101 ) 


( 2 . 102 ) 


{€ 0 /s)( V 

s 2 +2z/s+(z/ 2 +gu 2 ) 


(2. 103) 




2 2 


S +2Z/S+(Z/ 2 +OJ 2 ) 


•) 


(2. 104) 


<Vs)(- 


s 2 +^s+0J^_j 


s+z/ 


(2.105) 


Since all of the parameters of Z(s) 


are real, it is readily apparent that the 
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* * 

real condition is satisfied, i.e OJ Z(s) = Z (s ) . The six eigenvalues of 


Z(s) s U(s) are, 


€ 0<*>U 


\z' € ° S + s+(I/±jOJ H ) 


( 2 . 106 ) 


3 ~° s+I/ 


X =€«* + 


( 2 . 107 ) 


^4,5,6 ^° S 


( 2 . 108 ) 


while the real parts of the eigenvalues are 


„ l\ 1 € 0 <jufi{cr+i/) . _ 

R * X 1 J=€ 0 (T+ 2-5-^ 5- >0’ cr>0 

L i»2j 0 (a+z/) 2 +(aj+a) ) 2 ~ 

H 


( 2 . 109 ) 


R e[X 3 j = ^ + -^^->°, cr> 0 


(2 , 110) 


Re [ X 4.5,6 

Therefore, the cold magneto-ionic medium satisfies the positive real condition, 

2 10 Causality Property 

The energy condition is not sufficient for a constitutive matrix to repre- 
sent a realizable passive medium. Even though there exists a causality condition 
for circuits, it does not play as important a role as the causality condition 
does for distributed systems, the reason being that the velocity of propagation 
is assumed to be infinite in circuits,, In distributed systems in which 


= fi. cr>o , a>0 


(2.111) 
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characteristic lengths may be large compared to a wavelength, causality implies 
that the wave front can travel no faster than the velocity of light in a vacuum. 
Hence, causality for distributed systems should be stronger than causality 
for circuits. And one might even expect that if the velocity of light in a 
vacuum were mathematically forced to approach infinity the causality condition 
for distributed systems would approach the causality condition for circuits. 
Thus, the problem is to use the causality condition to derive extra necessary 
conditions upon the constitutive matrix for realizable passive media. 

Consider an infinite homogeneous passive medium in which a Dirac delta 
electric or magnetic current source both in space and time is placed at the 
origin of a coordinate system. Then causality states that for the field 3P (r^t) , 

CftTt)=0 for t ~ < 0 (2.H2) 

where c is the velocity of light in a vacuum. 0^ (r^t) may be represented as 
a Laplace transform in time, 

V i°° 

(2.113) 

where <J q is a finite constant for which F(r,s) has no singularities in the half 
plane Re s > CT q > 0. Evaluate Equation (2.113) for t - r/c < 0 by closing the 
contour in the right half plane. Then we have 

^{T.t ) = (H/ 2 irj) JPim /f( 7 ! s) e~c $(t c> as = O i t --=-<0 (2,114) 

R— *.0D''C 

where C R is a semicircular path in the clockwise direction of radius R. Assume 
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that the integrand satisfies Jordan’s Lemma, i e , 

|F(T,s)e | < when lsl>R 0 for Res>0£ 


( 2 , 115 ) 


where K is a positive constant and M(r) is some vector independent of s, Then 


nr,s)|<^-i^ ReS 


( 2 , 116 ) 


when jsj > R^, Re s > CT^, for all r Therefore, we expert that for any reali- 
zable medium the Laplace time transform of the field due to a Green's (Dirac 
delta) source should obey this necessary condition Note that it is easy to 
verify that a vacuum medium satisfies this condition, 

Since it is difficult to find a more explicit condition upon the consti- 
tutive relationship for an arbitrary medium, assume that the medium is isotropic 
Also assume that the elements of the constitutive matrix are ratios of two 
polynomials in the variables s, and k. Let, 


det(0-su)= D(s,k) 


(2 117 ) 


del y sqts,k)/D (s,is) 


(2 118 ) 


(0-3u)‘= -=^3I 

D(s ,k) 


( 2 , 119 ) 


where D^Sjk) and D 2 (s,k) are polynomials in the variables Now we have, 


(0 ~ s y )~'= Jk (s ’ k) ^ (s > T) 

M D 2 (s,k)D(s,k) 


(2. 120) 
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with the elements of D 2 (s,k) N(s,k) as polynomials in the variables. Therefore, 
the field due to an arbitrary Green's source (r^t) = C 6(r) 6(t) is 


^rrr D 2 (s » k) ^ s ^)c A -jir.r, 3 , 

FKn.umJfJ e dk <2121) 


€0 

:( 277 i 3 [D JsjV)N( Sl jV)]c/yT ^ 

L V© D 2 (s,k)D(s,k) 


since the elements of D (s,k) N(s,k) are polynomials. By integrating in polar 


coordinates over the two angles one obtains, 


F(r,s ) = \{27T) 2 D 2 (s,jV)N(s,jV) C 


Hf> kdk 


D 2 (s,k)Dls,k) 


( 2 , 122 ) 


And finally a contour integration yields, 


F(T,s) = |D 2 (s,jV) 


N(s, jV) C£l 


Residue 


D 2 (s,k)D(s,k)J 277-r 


-jk n (s)r 


(2.123) 


where k (s) is a root of D (s,k) D(s,k) = 0 that satisfies the radiation condi- 
n ^ 

tion Therefore, by comparing Equations (2.116) and (£.123) we find that for an 


isotropic medium, 


- 

Re jk„(s) > c'Re s , ls|>R 0 , Res > <T 0 , vn 


(2.124) 


where k (s) is a root of the determinant al equation 
n 
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det 


jofit)- s.y(s 7 k )j = 


0 


(2. 125 ; 


that satisfies the radiation condition. 

Now consider the special case of an isotropic medium in which both the 
permittivity and permeability are scalar functions of s only and of the form of 
the ratio of two polynomials. Then the causality condition implies, 

Re 


jk n (s) = Re s[^L (s)e(s)) / 2 J> c'Re 


S, isl >R 0 , Re 


s ><T 0 


■ 2 , 1261 


But, for both the permittivity and permeability, the positive real condition 
requires that the degree of the numerator minus the degree of the denominator 
polynomials is either, 0, -1, or -2 0 Therefore, for both the positive real 
condition and causality to be satisfied, the degree of the numerator equal a tne 
degree of the denominator for both jx( s) and £(s) r 

2 C 11 Example of a Bandpass Waveguide 

Finally, let us give an example of a medium that, if realized, can bo 
advantagenously used to produce a bandpass waveguide. Consider a lossless 
medium whose permittivity and permeability are, 

2 2 

* €(CU) = K — | % . M * Mo (2. 127 3 

cu 2 or 

Then s €(s) satisfies the positive real condition if and only if ^ is less 

than or equal to cj . But 

2 


( M/ IsI )' /2= <Mo k I' 2 2 c* 


(2 1 2 8 ) 
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implies that K > € qL Suppose K - It is obvious that the permeability 

satisfies the necessary conditions Now the medium inherently has a stop band 

between oo and co^ Fill a uniform waveguide with this medium. For a uniform 
i z 

waveguide the propagation constant Y is 


Y-- 


C Ll) 


gc 


^ H- o € . 


(UJ) 


(2 0 129) 


where to is the guide cutoff for a particular mode when the waveguide is 
gc 

filled with a vacuum. Then propagation exists for Y imaginary, or for frequen- 
cies oo such that 

UJ ‘ d fJL 0 € (u; ) - c 2 CU g C > o (2 . 1 30 ) 

2 . 2 2 *’ 2—12 

or ‘-CiJ — (jUgc ^ 0 c Figure 2 1 gives the regions of propagation. 

And the cutoff frequencies for the waveguide modes are 


OJL 




(Olp+ Cl)| c 






(2,131) 


Also observe that 


oo n 


and + are cluster points for the cutoff frequencies. 
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3. GENERAL FORMULATION OF THE SPECTRUM OF CHARACTER I STIC WAVES 


3 0 1 Field of an Arbitrary Source in a General Anisotropic (^,6) Lossless 

Medium as a Spectrum of Characteristic Waves 

Let us consider the problem in which the medium in Fourier transform space 
is described by matrix permittivity, €(k,cj), and permeability, y,(k,(jo), even 
though no known natural media has such a constitutive relationship. Also assume 
that the elements of the matrices are arbitrary functions of oj and k 0 With this 
formulation Maxwell’s equations display a great deal of symmetry 0 And the 
symmetry is capable of revealing much as far as the form of the expected results 
are concerned^ 

The Fourier transformed Maxwell's equations are 


- j K a E ( k . (jj : “ - $U)fi 0 N H ( k , UJi - J m ( K , QJ) 
-jkxHU.cji -- jcue 0 KE ( k ,uj) + J e ("k ,oj) 


(3 0 1) 


where r ~ ^ ^ and _ , ru : < 'jj ) Eliminating E in Equation (3 C 1) 

gives the equation for H c 


r, % 


A ~ K . 


v Oj) r 


~ !OJC 0 l^ m < k ■ 


(3.2) 


where 


M m ( k , LJ 





k x K J p ( k , 0 ) ) 


(3.3) 


Then Equation (3.2) may be rewritten as, 


G H (k,CJ)H(k ,Cl>) = -]CU€ 0 M m (k,GL)) 


(3.4) 
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in terms of the matrix operator, (^(kjU)), 

G h ( k , OJ ) = (- k x K 1 k x - N ) (3.5) 

The source-free solutions may be found from Equations (3.4) and (3.5) with 

the source term equal to zero. Alternately, one can find the eigenvalues and 
10 

eigenvectors that correspond to the problem 

-kxK 1 kxH; = X; N H; (i= 1,2,3) (3.6) 

2 

and then equate the eigenvalue to k . This means that the source-free equa- 

1 o 

tions will be satisfied for the propagation vector, k, on certain surfaces in 
Fourier space. These surfaces are prescribed by the equation det [G H (k,w)] = 0. 
In terms of the eigenvalues of Equation (3.6) the det [G H (k,co)J is 

(det N)(\ -k 2 )(\ -k 2 )(\ -k 2 ) or (det N)S,S 0 S„. 

10^0*30 1 Z j 

2 

(\.-k~) = s. = 0 is a portion of the dispersion surface. For a cold plasma 
(N = I , K ^ K(k)), S. = 0 is one sheet of the dispersion surface; however, for 
a warm plasma, (N = I , K = K(k)), = 0 (i = 1,2) may be more than one sheet. 

Assuming that the medium is lossless, i.e., K = K^, N = N^~, there exist 
certain orthogonality relationships between the eigenvectors IT. First, however, 
it is necessary to show that the eigenvalues are real. Define A as 
A = - k x K ^ kx. 

Lemma 3. 1 ; The eigenvalues of the matrix equation AH. = X.NH. (i = 1.2.3) are 

± i i > > 


real. 
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Proof z A H | - X j N H j 

A H j = X j H f N H j 
hf A t H [ = X* H 
Hj f A H ; - X* H^NHj 

(v A and N are Hermitian for real k and C 0 ) 

Eut Hj^AH; = X: ufNH; 

Therefore, (X; X^ ) I = 0 

! ) Xj = X? 

2) H j* ' N Hj = 0 

Theorem 3 t l z The eigenvectors of the matrix equation AH i = (i = 1,2,3) 

satisfy the orthogonality condition NH i = 0 for \ 

Proof: AH j = Xj MHj 

H: T AHj = Xj H^NHt 
Hj f AH j = Xj Hj'^N H ; 

(••• A,J\[ are Hermitian and Xj is real) 

But Hj^AH ; = X; Hj^N H; 

Therefore, (X; - Xj ) H j" N H; - 0 

Hence, J\l H j - 0 for X ; X j 

We observe that for lossless media, since the operator -k x K x N ^ is the 
- 1 — -!-► 

adjoint of -N k x K kx, then the set of eigenvectors, B^, are the reciprocal 
11 

basis to the set of eigenvectors, H , 


Hence , 
or 


37 


One of the eigenvalues, say \ 3 , is zero. This fact immediately follows 
from the relationship 

do t (-N kxK kx} = X| \ 2 X 3 (3*7) 

since the determinant of a product is equal to the product of the determinants 

of two matrices and since det (kx) = 0, The eigenvector, H^, corresponding to 

the eigenvalue X. = 0 then may be chosen as H = k. From this it is evident 

that H 0 is longitudinal and that E» corresponding to H_ is zero. 

3 3 3 

Choose the components of the eigenvectors, and to be polynomials in 
the transform variables £, rj , with no common factors. This is always pos- 

sible when the elements of K and N are the ratios of rational polynomials of 
the wavevector k and uj. Then the identity matrix, I, in terms of the eigen- 
vectors is 

3 

I = X ( H ?NHj r' NH i H f (3.8) 

i = I 
3 

I = Z l H ! M H i f 'Hi h|N (3.9) 

i = l 

The operator of Equation (3.8) operating on a vector splits the vector into its 
components that are parallel to B^. Similarly, the operator of Equation (3,9) 
splits the vector into its components parallel to IT. Equation (3.8) will be 
used in the interpretation of the field due to an arbitrary source 0 

Both (^(£, 0 )) and its inverse g” 1 ^,^) can be expressed in terms of the 
eigenvalues and eigenvectors of Equation (3.6). 
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G h ( k ,CJ) = Y, S; (h| NH| f' NHj H* N 


(3.10) 


G‘ H (k ,CJ) = X Sj'(H| N H | )"' H; Hf 


i = i 


(3.11) 


To verify that the expression for the inverse of G H (k, w ) given in Equation 

(3.11) is correct, multiply (^(k^G" 1 (k,u>) on the right by the vector NH . . 

J 

Then 

3 

G h G H NH j = G hZ s f 1 ( h | NH; )"‘ Hj h|nH- 

1=1 J 

= Sj 1 G h H j 

= Sj Sj NHj (v Gh Hj = Sj NHj ) 


= NHj 


Since any vector may be represented in terms of the three eigenvectors, 

H i (i = 1,2,3), Equation (3.11) is true in general. The solution for the 
magnetic intensity, H(r,a)), for the arbitrary source is obtained by the 
inverse Fourier transform of H(k,co) deduced from Equations (3.4) and (3 0 11). 

00 3 

H(7,aj) = ■■ )U)€ 0 [2w)~*fff f, S" i '(HfNH i r l H i H|M m (k,CU)e jk ' r d 3 k (3.12) 

00 1-1 


Multiplying Equation (3,12) by jx inside the integral sign gives the expression 
for the magnetic flux density, B(r, w ). 

k (3.13) 


00 


B(r,CU) = -jCU Mo € 0 (27T) 


-3 


3 

Is' 


(Hj + NH,- r'NHiH|M m (k,Ct))e Jk;r d 3 


- 00 
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Recall that Equation (3.8) operating on a vector M (k,u)), splits the vector 

m 7 7 

into its components that are parallel to the magnetic flux density B that 

3 ’ i* 

corresponds to the eigenvector, H.. Therefore, M m (k,OJ) = £ M mi (k,CJ) 

1 i=l 


M mi (k,CJ) = (H|NHj ) ' NHj h{ M m ( k , ( a ) ) 

^ 3 

B(k,aj) = X Bjtk.o;) 

i= I 

Bj (k,CJ) = (h{ NH j f'NHj h| B(k,6t>) 

Using these facts result in 

co ^ 

B(T,oj) = -iW 0 € o (27r)~ li [ff X S i , M ml (k,aj)e" jk r d 3 k (3.i4) 

JJJ j = | 

-oo 

Because of the orthogonality condition, Hj NHj = 0 (H|Bj = 0) Xj 7^ Xj , the 

component of B(k,o)) parallel to B is entirely due to the component of M (k cu) 

1 m 9 

parallel to B i , i e e. , 

Bj (k,a>) = -\LJfi 0 € 0 S]' M mj (k,CJ) (3.15) 

Denote the inverse Fourier transform of Sj 1 and (H* NHj) - 1 by Q. and G. , res- 
pectively. Then 

Bj ( r,0J) = -jCL»/i o e o Qj * M m j(?,aj) (3.16) 

and 

3 

B(f, 0 ))= 0 j * M mj ( 7 ,CU) 

i = l 


where 


Similarly 


where 


(3.17) 
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M . (r ,w) , the inverse transform of M . (k,cj) is 
mi ’ mi ’ 

M mi (?,0J) = /oGj * B| ( jV,CU) s * j V, CU)S * M m (T,CJ) (3.18) 
=>' 0 'Gi * Bj ( jV,0J) Hj f (jV,CU) M m (T, CJ) (3.19) 

The prime symbol is used to emphasize that the quantities are eigenvectors , 
not components of source fields. The reason for the choice of the normaliza- 
tion for and now becomes apparent in that the inverse transform of H'(k,co) 
and B^(k,oo) are H'(jV,oj) and B^(jV,co), respectively. Also, because of the 
isomorphism between polynomials in Fourier space and partial differential opera- 
tors operating on the Dirac delta function 6, an orthogonality condition equiva- 
lent to HjNHj r 0 1 Xj7 Z "Xj , exists in x,y,z space, i. e. , 

Hj + ( jV, GU) S * B[(jV,cu)S =0 for XjT^Xj (3.20) 

or 

Hj + ( jV,C(J)Bj' (jV,CU)8 - 0 (3.21) 

If in Equation (3.1) H was eliminated instead of E, analogous results would 
occur. However, invoking duality produces the same results in a more enlighten- 
ing manner. Duality implies that one can replace the quantity (I) by the quan- 
tity (II) in the previous formulation to give the desired results. 
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(I) 

E 
H 
J 

e 
J 

m 

e 

o 
N 

K 

Using duality, with the dual operator S , Equations (3.2) through (3.5) result 
in, 

(-kxN ! 1Tx-koK)E(ir,CU) = - )CO/i 0 M e (k t (x)) (3.23) 


(ID 


H 

-E 

J 

ra 

-J (3.22) 

e 

€ 


^o 

K 

N 


M e (K,o»=-«S = J e (k,aj)+(i/jaj Mo )(-jkx)N' l J m (iT,a;) 0.24) 

G E (fi,u;)E(K,aj) = -ja)/* 0 M e (t,aj) 0.25) 

G e ( k ,cj) = £) [g h (K, 0))]= -k*N*kx-kQK (3 26) 

Since the operators det and JB commute, det [G^,(k,cj)] is easily seen to be 

(det K) [S(X,-k*)][.B(X 2 -k*>][i>(X 3 -k*)] 

- 0; therefore, s [ 53 } 53 . At this point, we will prove some relation- 
ships between the eigenvalues and eigenvectors, that are necessary in order to 


continue the discussion of the electric field using- the duality principle. 
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Lemma 3.2 : The dual of eigenvalue X; is Xj ,i , j * 1 , 2 , i^j , if 

•£[xi]#Xi. (3 3,£»[xi]#x 

Proof: Assume that an eigenvalue, X^, and its eigenvector, IT, satisfy the 

equation 


G H (K,CU,Xj) H j * 0 


(3.27) 


The dual of Equation (3.27) is 

J B [G H lk,0),X )Hjj 8 G E ^k t CJ,^^X,jj Jb [hJ* 0 0.28) 

Therefore, £) [X. ] is an eigenvalue of G (k,w,\)E = 0 with^ [H ] as its 
* i £ i 

vector. Since 

Idet N)det[G E U,a),X)] =(det K) d€t[^G H (k,a),X) j = 0 (3.29) 

JE> [xl is also an eigenvalue of G H (k,w, X)H = 0. Nowjg) [X A ] £ 0 since it is 

assumed that X. t 0. Therefore, either £) [X. ] = X., j = 1,2. But, since 
X ^ J 

&[\] t £l\] = i,j = 1,2, i * j. Note that if £[\] = \ for 
either i equal to one or two, then the other nonzero eigenvalue must also be 
equal to its dual. 

The duality relationships among the eigenvalues and eigenvectors are 
summerized in Figure 3.1. 

Now with the aid of Lemma 3.2, we are able to discuss the implication of 
duality upon the dispersion surfaces. First let us note that the dispersion 
surface for the electric field is not necessarily equal to the dispersion sur- 
face for the magnetic field. By this is meant that the zeros of det Gg are not 
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XjsO — GjEk,0» f X)— H 3 =k M.E. ► E 3 =0 Magnetostatic 


Electrostatic H 3 = 0 M.E. E 3 s k G E $/n l X) — X 3 = 0 


Figure 3.1. Duality relationship among the eigenvalues and 
eigenveetors of_^Gg (k,w,^) and Gg(k,w,\). 
(GgtkjOj,)^) = - k x N -1 kx - XK; (^(kjWjX.) = 

-k x K -1 k x - XN;^ = dual; M.E. = Maxwell's 
Equations) . 
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necessarily equal to the zeros of det Gjj. A warm plasma is an appropriate 
example. Hence the phrase, the dispersion surface of an electromagnetic field," 
is not precise. For the dispersion surfaces to be equal would imply that the 
zeros of det € equal the zeros of det jx. In truth, the relation between the two 
dispersion surfaces is duality. Lemma 3.2 can be used to explain a finer struc- 
ture among certain sheet(s) of either dispersion surface. Consider the disper- 

3 

sion surface for the electric field, i.e,, det G_ = (det K) 'll S. = 0. The sheet(s) 

E i=l 1 

defined by S 3 = 0 is equal to the dual of itself. If the sheet (s) S. = 0, 

(i = 1,2) is not equal to the dual of itself, then the dual of sheet (s) S. = 0 
is equal to the sheet (s) S = 0, (i 4 j, j = 1,2). This characterist ic is exem- 
plified by the general time-dispersive uniaxial problem of Section 4. Similar 
results apply to the dispersion surface for the magnetic field, det G = 0. 

After this discussion of the implications of duality upon the dispersion 
surfaces, eigenvalues and eigenvectors, the procedure to obtain the electric 
field from the first part of this section by duality should be self-evident. 

From the previous formulation, the field ^=^6^ could be found in 
three ways: (1) by the equation G u (k,cj)H(lT,co) = - j<o€ M (TT.co) and Amperes law. 

(2) by the equation G^ (IT,cj)E (ET,u) = - jwp- o M e (k,o>) and Faraday's law, and (3) by 
the equations G H (lT,u)H(k,w) = - jw€ Q and (£,w) = - jwp o (lT,co) . 

Only the last method is symmetric in the field components. All three methods 
depend on either G^ or G^,. That is, in all three methods. Maxwell's equations, 
with an assumed form of constitutive relations (D = €E, B = jpl) , were reduced 
to an equation envoi ving only one field component, E(lT,oo) or H(£,o>). However, 
this relatively simple reduction may not always be possible. If the flux 
quantities, D and B, are linearly related to both the field intensities, E and 
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H, then the simple reduction is not possible. An example of such a coupled 
constitutive relation occurs for the fields in a moving medium. In general , 
the linear constitutive relationship may be written as 


F f = UF where F = [D,B] and F = [E ,H] . 

Thus, it seems desirable to formulate a method to find the field due to a source 
in terms of the characteristic fields in a symmetric manner that envolves both 
E and H, and that circumvents the difficulties that arose in the previous formu- 
lation. 

3.2 Definitions of Symbols 

Let us define some symbols of quantities (vectors, matrices, functions) in 
the Fourier transform domain. 


F = 


Y 

H 


(six-vector of the electromagnetic field elements) 


(3.30) 


C = 


(six-vector of the electromagnetic source elements) (3.31) 


ml 


ro -j&i 

0 '[jiTx o J 


(3.32) 


U 


(sixth-order constitutive matrix) 


(3.33) 


7 Tj. dT,u>) = 0 - jcoU 


(3.34) 
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F f = 


B 


= UF (six-vector of the electromagnetic flux elements) (3. 35) 


I (sixth-order identity matrix) 


(3 36) 


F i ,v i ^ the ei ^ envector s and eigenvalues corresponding 

to the matrix equation OF = v UF ) 

i i = i 


(3. 37; 


F f i = M F i 


(3 38, 


f i“ -I + 

ifMF-UFj) UF ( F( 


{3 39 


If a (I i , ) f =(F}yF.j H F.Ftu 


(3. 4G 


V fi “ I i V an arbitrary vector; V component of V 


parallel to F ) 


(3 \ 


V i = I i v <- v an arbitrary vector; V i} component of V 


parallel to F ) 


(3 42 


3,3 Characteristic Waves 

By definition, characteristic waves are the fields which are described t 
the eigenvector, F^ , and its associated eigenvalue, v.^, that are the solut • 
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the equation 


0F 7 




(3.43) 


Since 0 and U are sixth-order matrices, there will be six characteristic waves, 
i.e. , the index i will run from one to six. When the eigenvalue, is equal 

to jcx>, then the eigenvector, F^ , will be a source-free solution to Maxwell’s 
equations. However, the field, F^, corresponding to v = v, may not encompass 
all of the solutions to the characteristic equation. Field vectors for the 
characteristic equation exist, if and only if 


det 



(3.44) 


But 



6 


(dety)II(l/j- y) 

i=i 


(3.45) 


Hence, other solutions will exist when det U = 0. It is not necessary that 
the determinant of U equal zero identically. In general, U will be a function 
of the transform variables u> and k. Also, U will be dependent upon other 
parameters. Thus for certain oa, k and other parameters of U, the condition 
det U = 0 may be satisfied. 

At this point it should be emphasized that the eigenvalue problem that 
concerns us is in distinct contrast to the usual eigenvalue problem. The usual 
eigenvalue problem deals with the question: Given a medium, normally isotropic 

and homogeneous, contained in certain boundaries with "walls" that may in gen- 
eral be described by an impedance, what are the source-free solutions or modes 
that may exist? These eigenvalue problems encompass both the bound case 
(cavities) and the unbounded case (waveguides). They are in the time-space 
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domaixr The modes then are a superposition of the source— free waves that may 
exist in the unbounded medium in such a way that the boundary conditions are 
satisfied For example, in a vacuum filled perfectly conducting rectangular 
waveguide, four appropriately chosen plane waves comprise a mode. In a per- 
fectly conducting circular waveguide, an infinitely nondenumerable number of 
appropriately chosen plane wave comprise a mode, The eigenvalues in this type 
of problem are largely determined by the boundaries. In our problem, however , 
the medium is of infinite extent. Also the domain of consideration is the 
Fourier domain in contrast to the real s p ace * t ime domain. The eigenvalues and 
eigenvectors are entirely a function of the matrix operators 0 and U. The 
operator 0 is a result of the fcrrr of Maxwell’s equation only and is indepen- 
dent of the medium The operator 0, however, is the constitutive relationship 
and hence comprises the entire electromagnet ic description of the medium. 

Thus, it can be said that eigenvalues and eigenvectors are a function of the 
medium only. In fact, the constitutive relationship U can be completely ex- 
pressed in terms of the eigenvalues and eigenvectors. 


Proof 


1L= ( Tied) 


MJ 


(3. 46) 


whe re 


U = 


o+jcju 

SfTufT'uf. fIu 
. = 1 ' = ■ « = 



(3.47) 


(3 r 48) 


(3,49) 
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Since for i = 3 and 6 both v i and OF^ are zero, the limit must be taken for 

these terms. It must be emphasized, however, that knowledge of and F. at 

one frequency is not sufficient ; and F^ must be known as functions of ^ 

and co. Equivalently, a knowledge of the characteristic waves at every point 

on the dispersion surface in four-space (k,oo) completely determines the 

electromagnetic properties of the medium. 

Certain properties of the eigenvalues and eigenvectors may be determined 

by only knowing the symmetry properties of the operators 0 and U. It is easily 

shown that the operator 0 is skew-Hermit i an for real k, i.e., 0 = - 0^ Assume 

that the medium is lossless. This is equivalent to saying that the operator U 

is Hermitian, i.e., U = . With these two facts it is easily shown that v. 

is imaginary. This follows from the fact that (fT OF ) is imaginary (F^ UF ) 

i i *i=i 

is real and that = (F^ 0F^)(F^" UF^) \ Also, there exists an orthogonality 

property among the eigenvectors. The conjugate transpose of the equation 

0F i = v i F t is the e Q uation OF. = v.FtuF.. Use has been made of 

JiJ J— x 

the symmetry of 0 and U, and that is imaginary. But we independently know 

that F+ OF = V.FTUF.. Hence, (v . - v.)F*^UF.. Therefore, for v ¥ v * 
ji ij=i l J J- i 9 i j * 

FtllF. = 0 or F^F _ . = 0. 

J 1 J fi 

To better understand these characteristic fields for an arbitrary Hermitian 
constitutive relationship, it may be profitable to digress for the moment to 
correlate the characteristic fields from the six-vector and previously considered 
three-vector methods. Naturally such a correlation implies that the constitu- 
tive relationship is "diagonal," i.e., U = 

decomposition into three— vectors possible. 


| 0 
0 y. 


since only then is the simple 



3.4 Relationship Between 77} (v) , GgCX), and GhOO 
It has already been established that 

G E ( X ) = “ k x Kx - X £ 0 £ 
G H (X) = -kx€ 0 € 'k x-\ fji^ fx 

and 


Further define 


2fU/) 


ve -jkx 

♦ 

for u = 

I 0 

1 

X 

1 

if* 


0 £ 


G e (\) z fj-Q G e (X) 
Gh(X) = £o'g e (X) 


Using these definitions it may be shown that 




G E‘-Mo e o l/2 > 0 

0 ^-ho ( o l/Z) 


Since 
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and 


g e 0 



0 


( 

xX 

o 

o 

_ 1 


0 

I 



i o 
O G' h 


(3.57) 


then. 


Ge o 
o g' h 


det /^(-i/) = det ond det 

Therefore, 

det ^(-l/)u'^(l/)] = (det u’S^det 


= (det G' E )(det G„) 


(3.58) 


= det G e (-c 2 I/ 2 ) det G' h {-c 2 v 2 ) 


Or 

det 2/) - ^ j(det det GeI-c 2 v 2 )j |(det €) det Gh( - c 2 V ijj 1 (3.59) 

This expression can be simplified further by using the following lemma. 


Lemma 3.3: (<j e f ^)d e t ) = (det |)det G^{~C Z J/ 2 ) 


Proof: Define 


// . 2 2. _ ~ I . 2 2 

G e (-C V ) - A/1 A + c V € 


(3.60) 


for an arbitrary matrix A 


r hen 


^ 2 

A £'g"(-C 2 i/) = gVc (3.62) 

det Adetf detG^ = det G^detyLX 1 de) A (3. 63) 


Or 


-2 2 


det^detGjft-C Z/ 


) = det£ detG^t 

n 


C Z l/ Z ) 


( 3 . 64 ) 



( 3 . 65 ) 


Now the determinant of 7?^(v) can be expressed as 

det 7Tj[l/) - + detyy, detG^f-C^Z/ 2 ) : + det€det G^(*C^i/^) (3.66) 

I' he ^ign ambiguity occurs since the square root has been taken. However, since 
tin determinant of 7T^ ( v ) is set equal to zero to find the eigenvalues, the sign 
iiiuient is of no great importance at this point. 

A few observations concerning the preceding mathematics are in order, par- 
t icularly about the relationships between the eigenvalues and eigenvector of 


the three-vector and six-vector methods. One apparent discrepency is contained 
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in the question: Why does the three-vector and six-vector methods result in 

three and six eigenvalues and eigenvectors, respectively, particularly in 

light of the fact that both methods represent Maxwell’s equations for the same 

medium and both sets of resulting eigenvalues and eigenvectors are sufficient 

to completely represent electromagnetic propagation in the medium? The answer 

to the question is found in Equations (3.55) and (3.66). Both Equations (3.55) 

and (3.66) indicate the relationship between the eigenvalues of the two methods, 
-2 2 

i.e., \ - - c v . This shows that corresponding to each of the three eigen- 
values X there exists two eigenvalues v = jc(\) 1/2 . if f = [E,H] is an 

eigenvector, i.e,, 7)T[ (v)F = 0, then from Equation (3.55), E and H are eigen- 

vectors of G^, and G^, respectively. Conversely, a pair of eigenvectors E and 
H, corresponding to the eigenvalue X of G^, and G^, is also an eigenvector of 
*1 (v) - Moreover, the six eigenvectors are not completely unrelated as one 
might surmise from the fact that corresponding to every eigenvalue v., there 
exist another one, There is in fact a relationship between eigenvectors 

corresponding to v and -y .. 


Assume that corresponding to the eigenvalue y, F is the characteristic 
field, i.e. , 


OF = v JJF 



(3.67) 


y 

Now we claim that the characteristic field, F, corresponding to the eigenvalue 

V 

-v is F = [E , -H] , i.e., 

v v 

OF = - v UF 


(3.68) 
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This is easily seen by distributing the negative sign of -H with the preceeding 
matrix elements and the multiplying on the left by , The result is 

Equation (3.67) which was assumed to be valid. 

The following table will summarize the relationship between the character- 
istic fields of the two methods. 


TABLE I 


X , E, H, 


X 2 E 2 H 2 


F l = 


E, , H| 


‘V\ F i = 


E , H 




” “ 


V 

- 

v 2 

f 2 = 

E 2 ,h 2 

-v 2 

F 2 ' 

F 2 i _ H 2 


V * = 0 F t 


x, = o 


Ey-k H 3 --0 
E 3 r 0 M ” k 


E 3 =k,0 


y 3=0 f 6 = 


0 ,H 6 =k 


Three-Vector Characteristic Fields 


Six-Vector Characteristic Fields 


3.5 Completeness of Characteristic Fields 

Before any attempt is made to represent a field due to a source distribution 
in terms of the characteristic fields s it first must be established whether such 
a set of fields can represent the desired field, i„e 0 , does the characteristic 
fields form a basis for the vector space. Let us consider this question in terms 
of projections. 

2 

A projection 1 ^ on a Hilbert space, H, is simply an idempotent (1^ = 1^) 
linear transformation of H into itself. The range and null space of 1^ are 


M. = 

1 


I , X ■ X Z H j and N. = 


I , X = 0 : 


:x£h , 


respectively. Thus, the pro- 


jections I. partition the Hilbert range space. From this we see that defined 
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as <F+ UF i > 1 U (i = 1,2, ... 6) are indeed projections as defined above. 

Furthermore, they are orthogonal projections = 0 (i 4 j). In our case, 

the projection I partitions the range space into the space of vectors parallel, 

M^, and orthogonal, N^, to the characteristic vector F . If an arbitrary vector 

F is to be expressed as the sum of characteristic vectors, then it is obvious 

6 

that the sum of the ranges of the projections must span the space, £ M. = H. 

i=l 1 

Now the sum of orthogonal projections is a projection. Thus, it is not sur- 
prising that the sum of the projections is the identity projection whose range 

6 

is the Hilbert space, i.e., £ I. = I. Indeed, this is an equivalent statement 

i— 1 1 

as to the completeness of the characteristic vectors F^. 

Analogous reasoning applies to the set of projections |l f : i = 1,2, ... 6 1 

f f 1 

The main difference is that the set 1^ : i = 1,2, ... 6 partitions the 
Hilbert space into sets of vectors parallel to the characteristic vectors F fi 
not F^. Therefore, the latter set would be more profitably used with flux 
field vectors, whereas the former should be used with the field intensity vectors 

3.6 Spectral Re presenta tion of F i elds Due t o a Sou rce 

One of the most important purposes of representing the field due to a 
source as a spectrum of characteristic waves is to gain physical insight into 
the process of propagation in the medium. Most often the source problem is 
attacked in the Fourier domain by inverting a complicated coefficient matrix 
followed by a Fourier transform inversion. Since all of the sheets of the dis- 
persion surface will occur as singularities of the Fourier inversion integral, 
the integration process is extremely difficult. Also, the mathematics reveals 
little physical insight other than the singularities, i.e. , the dispersion sur- 
face which implies the source free solutions, are a prime determining factor in 


the source problem. 
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A number of questions may be asked concerning the relationships between 
the source free problem and the source problem c A very significant question 
is: what is the explicit relationship of each sheet of the dispersion surface 

and its associated characteristic wave to the source field? Also, in what 

manner does the source excite the spectrum of characteristic waves to form 

the source tic Id? These and other questions will be answered in the process 
of finding to spectral representation of fields due to a source. 

There is at least two ways of deriving the spectral representation of 

the source fields. The first and perhaps the simplest is to express the 

source in terns of fcnd assume that the source field may be expressed in terms 

of the characteristic fields. This is possible since it has already been 

established ; hat the characteristic fields are complete and may be used as 

a basis. Alsa, an important factor is the orthogonality of the eigenvectors. 

Actually the source will be expressed in terms of the characteristic flux 

fields F since the source itself is a flux as can be seen from Maxwell's 
5 i'P 

equations. That is, let 


6 6 

c ‘ I c fi = I F fi 

i= 1 


(3. 69) 


C is the c-mponent of the source that is parallel to the cnaracter ist ic flux 


field F . . 

1 l 


hen use the orthogonality condition to determine the coefficients 

Hence , 


6 . Therefore, F* C = 0. F f F . or (3. = (F f F_.) 1 F* C 

*1 5 1 r i 1 fl *1 1 fl 1 


C’ld^Fnf'F.^C 


(3.70) 


i = I 


S imi 1 ar : y , 


F ( k ,CtJ) = ^ Fj ( k ,CJ ) - 


i = I 


X a i F i 


i=l 


(3.71) 
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F^Ckjoa) is the component of the source field that is parallel to the charac- 
teristic field F^. 

6 6 

WT.a; IvXtfi F fi (3.72) 

1 1=1 i = I 

But 

/^Fj =(Z/.-jCJ)F fi (3.73) 

and 

F J + ^Fj = -jOJ)F i t F f; Sjj (3.74) 

where 6.. is the Kronecker delta. Therefore. 

ij * 

Q: (Z/| -jO)) Fj + F f j = F; + F fj (3.75) 

and 

( 2 ; =(Z/ i -jO)f (3.76) 

Equation (3.76) shows an interesting fact. The component of the source field 
parallel to the characteristic field is completely and entirely due to the 
component of the source that is parallel to the flux field F^. Moreover, this 
component is only affected by its own sheet of the dispersion surface and not 
by the sheets of the other characteristic fields. Now the component of the 
source field F^(k,co) parallel to the characteristic field F^ may be expressed 
in its spectral form, 


F i (k ,OJ) = ( i/. -jQjf '(FtFfjf V; Fj + C 


(3.77) 
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The total source field as a spectral representation becomes, 


Ftk.UJ) 4>ri“>'< F i fF «>Vi t c 

i=i 


( 3 . 78 ) 


The second method of deriving the spectral representation in contrast to 
the first is to directly express the matrix operators and 77| * in terms of 
the characteristic values and vectors without making any assumptions about the 
source or the field. If the first method had been determined at the outset, 
then it would give a motivation for the expression for //£ . Let us postulate 


the expressions for and in terms of the characteristic values and fields, 

and then show that they are indeed valid. The appropriate expressions for 
and 7f) are , 




( 3 . 79 ) 




( 3 . 80 ) 


The notations for I ^ and 1^ have previously been defined, 
sion for miuj) , it must be shown that 


To verify the expres- 


^(ja))v = ujj(z/-jca)i f jv 


for an arbitrary vector V. Since, however, the characteristic fields are com- 
plete it is only necessary to show that the equation is satisfied for vector 
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V to be any characteristic field. Assume that 

TT^ ijoilFj = u[X (i/, -jwiii] F ( 

L i= I -* 

Then 


( 3 . 82 ) 


7f[ ( jOJ) Fj = U(i/j -jcu)Fj ('•' orthogonolity ) ( 3 . 83 ) 

= (Z/j-jCt))yFj ( 3 . 84 ) 

= ♦ (z/j-ja))yFj (*•* ^ , (z/ j )F j = o) (3.85) 

= /^(jCL»)Fj (*•* ^(1/j) = O-J'jU) (3.86) 


Hence, the expression for^(jCJ)is verified. Perhaps the simplest way of veri- 
fying the expression for 7T[ '( ]CJ)is to directly show that^{j(ij)^ '(j<A») = I,i.e. , 

6 


1?l (ja))y -, [Z C2^ ( - ~ jCU) "*lj =1 


( 3 . 87 ) 


Now 


(jOJ){u H [ Z l I* J } = Z < -j^)’ l ^(jOJ)u" l I i f ( 3 . 88 ) 

6 

s Z tZ'j -jtjf'^tjcu) IjU ” 1 ( 3 . 89 ) 

1-1 ( V definition of I { and I-) 

6 

= Z u I ,• u _l (v P2{\(jj)i. = ua -jcdJuij ) 

r=) 


( 3 . 90 ) 
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(3.91) 


= I 


(3.92) 


Therefore, the expression for Tf\ ' ( jCd) i s also verified. 

To derive the spectral representation for the source field F(k, w), it is 
only necessary to use the spectral representation of the inverse of the field 
operator ( jCi))which was given above. Thus 


F(k,o;) = u " 1 X (z/j - jcu )" 1 i f f C 


(3.93) 


1=1 


= u~'X (z/j -jojf'Cfi 


(3.94) 


or 


F f (k,oj) = £ (^i-j^r'c fj 

i = I 


(3.95) 


Not only is this expression compact, but it is also quite revealing in physical 

insight. However, before the physical interpretation is given, let us show 

► 

one moi*e fact. Resolve the source flux field F^(k,w) into its components paral- 
lel to the characteristic flux fields. 

By definition the component of F (k,u>) that is parallel to F is 


.f „ .T 


fj 


F f j(k,w) = Ij F f (k,w). Using Equation (3.95), this yields 


F (k,co) = E (v - C . 

J i=1 x j 
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Now because of orthogonality - 6^j Therefore, 


F f j(k , oj) = (l/j C fj 


(3.96) 


and 


F f (k,o>) = Z F ti 

i = l 


(3.97) 


Again we arrive at the conclusion that the componehtT of the source field 

that is parallel to the characteristic flux field* F^_ , is solely excited by 

the component of the source that is parallel to the characteristic flux field, 

F . Now using Equations (3.69), (3.96), and (3.97), a physical interpretation 
f j V- ' - ' \ 

to the source field may be given. Simply stated, the Fourier transform of the 

flux field F (k,cj) is vi tiie ifum of the components of the source parallel to the 
f \ 

characteristic fields, C^., each divided by iis respective sheet(s) factor-of 

\ '%. sf**; 

the dispersion surface, S = (v.' _ jco). This fact is illustrated “in Figure 3.2. 
It is also evident that there exist -source distributions such that certain 
sheets of the dispersion surface and their respective characteristic fields 
play absolutely no part in the total field, namely, those source distributions 
which are orthogonal to the said characteristic f ields, i.^e. , C f s such that 
F^C = 0 for some index j. 

j 


In the real space-time domain, the flux field may be given as 

6 

: = Z * > r. r " 


(3.98) 


i= I 


where ^ f ;(7,t) = Q. * t? f j and Q.^ and£? f| are the inverse Fourier trans- 

forms of S -1 = ( v - 1w) _1 and C„., respectively. The convolution is a four- 
i i ±i 
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dimensional space-time convolution. In the space-time domain the source flux 
fidld ^(^t) in terms of the "characteristic sources" is illustrated in 

Figure 3.3. 

/*" } 

At fM^i^int the^ matrix Green* s function T is easily determined. By 
definition the matrix^Hreen* s function F (r, t) is an operator such that 


(3.99) 


inverse 


\ " f ; v '■ * 

Now the "characteristic source" is tv ?■:* « where ^ is the i 

\ f 

Fourier transform of the operator I = . Therefore, from Equation (3.98) and 

?? ' ^ ' •». 1 

• . . - *■ 

since the convolution operations are associative, the source ^lux field is 

* ’■> 1 -. 

\ 


\ ’ ' ■« 


( 3 . 100 ) 


Thus, the matrix Green 1 s funct^ioji in terms of bi|e characteristic fields and 

‘A ^ ‘ ■ \ 

% 

the influence of the dispersion surface (the operator is easily expressed 

1 \ \ 

in terms of the characteristic fields) is ‘‘V v 


i=l 


- (3.101) 

/ 

/ ■ H *. jS 

H- A • - ' 


Alternately, the matrix Green’s function could have immediately beenxdeterAiined 

^ -If 

from the inverse Fourier transform of (jCJ) which is equal to E 1^. 

It should be emphasized' that these results are general and that the only 

• ■ I , ' ; : • • 'r- ’ C **.i ■ ' '■ ‘ r ' '• . > • , • . •; . ij - 

conditions that have been made a^re that the medium be homogeneous, linear and 


lossless, i.e., have a Hermitian constitutive relationship* 


65 


3 . 7 Bquilvalence Between Formulations 


Again it is instructive to show the equivalence between the six- 

f- 0 ** 

and three-vector methods for the source problem. Let U = 

° fct 

U = U. The off "diagonals" of U must be zero for otherwise the~s , imp] 
coupling between the E and H vectors does not apply. 


n 


F = C 


Multiply both sides by U * . 
But 


-7il u"c = -jojc' = -jcL> 



Je_ 


v 

jOJ 

_ 1 - 1-» 

= -jO) 



jjj kx€ kx 

Mm 


Therefore, 


= -joic* 


F (k,0)) = 




! - 1 
6 


3f_"' = X t Z^i-jCU)” 1 X|U _I 


1 = 1 


Therefore, 

66 

F(k,a» = - X -jwf , (i/ n +jajr , iiy“ l ui n y‘ , H^C / ) 

i*l n=l ~ ~ 


vector 
and 
e de- 

( 3 . 102 ) 


( 3 . 103 ) 


( 3 . 104 ) 

( 3 . 105 ) 

( 3 . 106 ) 

( 3 . 107 ) 


( 3 . 108 ) 
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= ■£ U/f + cuVljiA-jidC') 

i = l 


(orthogonality) (3.109) 


© 

F(ir,tu) = -ja) Mo e 0 I (X, -k„r li uc' ,(••• Xp -v, ho e 0 ) w.n m 


Since \ = \ 3 \ = X. and C' is independent of an index, we may sum these 

14 2 5 

terms. Therefore* 


FjU,aj) + Fi +3 (k,CL)) = -jCU Mo € 0 (Xi-kof , (Ii + Ii + 3 )y' , C / ,(i-l,2) (3. ill) 


ii = (Fi + yFj f'FjFty 


(3. 112) 


(Fj UFj ) = (Fi + 3 yF i+3 ) 


(3.113) 


for i = 1,2 since F. _ = F. 

’ 1+3 i 


Therefore , 




I, + li + 3 = 2 (F, UF, ) o t y 


E,E| 0 


(3.114) 
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And 


c. 

F(k,oj) = -2jcL» Mo € 0 Xu i -kor , (Fi + yF i r l 


1 = 1 


E t Ej 0 


0 HjH; 


C' 


€ o 1 (^-kJr'lF+yFjf'F.F^C' 

i = 3,6 


(3.115) 


For i = 1,2 


For i = 3 F, 


( Fj + UFj ) = 2{E *€Ei ) = 2(HjVH; ) 


= [e 3 .°] 


and for i = 6. Ft 


= [o.Hg] 


Using these simplifications, Equation (3.115) reduces to that of the earlier 
method through and G^. Hence, 


E(k,OJ) = -j<U/i 0 € 0 X(X,-kJf^Ef CEj’f'E, E, + 

i = l 


M, 


(3.116) 


And 


H(ir,a» = -ja>Mo«o I U i -k*r'(H i VH i r l H i H> m (3.im 

1 = 1,2, 6 

3.8 Equ ivalence Between Formulations By Direct Addition of Terms 

An alternate and perhaps more illuminating correspondence between the two 
formulations may be made by the direct pairwise addition of the terms corres- 
ponding to and -v^. 
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F j(k,OJ) + F i+3 (k,cu) = (l/j -joj) ' Ij y H C + (-Z/j -jCxjf'lj u"'c (i=l,2) (3.na) 


= (^i -jcuf '(F^UFjf ' FjF^C - (^i + jCU) l (Fi t yF i r l FiFi + C 0.119) 


: i + y Fjf 1 U'j - jOJjVj Fj f -U/,+j(u")F i Fi t C 


( 3 . 1 20 ) 


(Fj UF; 


EiE 0 


[V f+OJ 2 ) _ 0 HiHit _ (i^f+CJ 2 ) _ H * E ^ 0 


o EjHj 


C ( 3 . 121 ) 


mm a 3 , 4 : 


^ i j CiJ El ] (jj ^ 


-jkxEj = -Z/j/xHj 
Z/j H| = j^'kxHj 

V\ Hj + = “jEfkx^ 


I/ i H i 1 'j m = iwE i t (- z Lkx g ' l J m ) 


J/iElje = jOJHt ( 


- j k x H j = l/jf E, 

Z/j E i = -j|“'kxHj 
Z/, E* = jHi + kx € -1 

= j oj H i f ( Jjkxe''j e ) 


Proof : 
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Using the two lemmas it can be shown that 


E i E i t (-s;kx g "j m ) 
H, Hj f ( Jjkx|''j e ) 



Ej Ej f 0 



'j m ) 

= 2 jo) 

0 H,H, + 


_( <b*«r 



2 V\ 

0 Ej Hj 

f 

C = 2 jOJ 


Hj E ; 0 



Placing this in Equation (3.121) and collecting common factors we get 


Fj(k,a;) + F i+ 3 (k,cu) = (FfuFj )"' -I jaj - 

= (1/f+aj ) 


EjEj + 0 

0 H j Hj f 


C + 


J m) 

( ikx|"j e ) 


(3. 


(i = 1,2) 


However } the term in the braces is just C* . 


Therefore y 

Fi(k,aj)+Fi +3 (k,aj)=-2jai /io 6 0 (x r kor l (Fi + yF i r l 


EjEi + 0 
0 HjHj + 


C' 


(3. 


(i* 1,2) 


Since X. = 0 and H = 0 

*5 o 


F 3 (k,CU) = (Z/j -jeuf'tFg.UFjf ' f 3 F 3 + C 


. 122 ) 


123) 


124) 


125) 


(3.126) 
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-jOJ/i 0 € o( X 3- k O> Wy p 3> ' F 3 F 3 C 


(3. 127) 


F s <k,a,l*-i^.e„(X 3 -K 2 0 )- , ( F:yF 3 r'F 3 F 3 t C' <»•«•> 


Similarly , 


F.(k,CJ) = 


-jGU)i 0 F 0 ( ^6~^o*^ F 6 yF 6 ) _l F 6 Fg C 


(3.129) 


= -icO / x 0 eo()'6-''o»"( F 6 UF.r'V«C' , (v He II k> (3.130) 


adding the terms gives the same equation which was found by the previous 


method , i. e. , 


Flk,CU) = -2ja)/i 0 € 0 Z Ur k o> V^UF,) 1 0 H H t c ' 


EiEj 0 


(3.131) 


-iw M0 « 0 I Ur* 0 > (F, yF|> f i f i c ' 

i = 3,6 


which reduces to 


E Ik, CD) = -jOL)/x 0 € 0 X U| ' lE i + I E '< ) E i E i f|V, e 


(3. 132) 


H(k,aJ) = -jCJ Mo ^o Z lX i -kJ)" , lH, + ft H i )" , H i HfM m < 3 - 133) 

i= 1,2,6 
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4. SPECTRUM OF CHARACTERISTIC WAVES IN A GENERAL 
TIME-DISPERSIVE UNIAXIAL MEDIUM 


The general time-dispersive uhiaxial problem is described by a matrix per- 


" >,T.. - i 'X ■■ - ■ - . fN, 

0 

0 ■ 

Kj 

0 

0 " 

meability and permittivity oi the form ZX = iX 0 

N, 

0 

and 0 

K| 

0 

[o 

0 

No 

[0 

0 

Kc 


respectively. Maxwell T s equations for a lossless medium are 


Curl ECD =-jaiyLZ. 0 NHCD - J m (r ) 

Curl H(7) = jCJ6 0 KE(T) + J e (T) 

(4.1) 

Eliminating the electric field gives the equation 



(4.5) 
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Explicitly G (k) is in terms of the wave vector k - 


_ ; l o £? 2 ~ l '££ 

g h oo = _L oC 7 ? ( L «C +L -C _k o N .> -L77C 

‘ L '££ -L,^C H-iC 2+ L|^-4n 0 ) 


(4,6) 


The characteristic waves are obtained by determining the solution to Equation 

(4,4) when the source M (k) is zero. This implies that the determinant of 

m ■ 

— ♦ — £■ 

G (k) is zero for nonzero solutions. Denote the determinant of G (k) by D(k), 
H H 


* - *0 1 L o T ./° 2+ L ,T,£' 2 - Ko > ( L , T 0 p* - L ,T, £ 1 ko 1 I N ! 


(4.7) 


where 


2 2 I 

e yO = (£ + Tf )* 


moxe , and by 


S,= (LoT.^+U.T, ^ -k 0 ). 


S l =U,T o/ > 2 +L,T,£ 2 -k=), V k ' 


(4.8) 


where T* N 


and det N s |N| 


-* 4>* 

A*-kxKkx 


73 


AH, * XjN Hj (i = l,2,3) 


<4. 10) 


or 


( A— X, N) Hj = O 


(4.11) 


We claia that X.. is real and H. NH. =0 
i J i ’ 


Proof : 


h|ah,*X,h|nhj 


(4.12) 


(HtAH,) + = X, (H^NHjj^ 


(4.13) 


+ V* + 

HjAH, = A, H, NH 


since A and N are Hermitian (4.14) 


(X,-X|)h|nh, 5 0 implies Xj real if H,NHj^O 

H|AHj«XjH|NHj 

h|aHj s XjH^NHj 

hJaHj = XjHj^NHj 


(4.15) 


(4.16) 


(4.17) 


since X_ real and A, N Hermitian (4.18) 


<Xj~ Xj)(»-i^N Hj) = o implies H^NH, = 0 if Xj^Xj 


(4.19) 




X 


Figure 4.1. One quadrant of the three-space dispersion 
surface for a general time dispersive 
uniaxial medium. 


75 


The eigenvectors are. 


Hj = kxz, H 2 =N l tx'kxz, Hyt 


(4.20) 


Then 



G „(lT)=X S i( h T n h { )' NH f h|n 

i=l 


(4.21) 


(4.22) 


^(Tl^sV.HjNH-, 1 H ( . hT <4.23, 

H+NH 1 .| V f HtNH 2 =N>-;(N |/0 S +N 0 CV 2 . H+NH 3 =(N, p\ H g g) (4.24) 


h,h| MTTxz^dTxz) 


T? 2 "^ 0 

C 7 ? £* 0 
0 0 0 


(4.25) 


HjHg - (n'ITx kxJiWkiliuf) 


"X 


<?C 

VT -N 0 N t 7jp% 

-n^vpr 4? 


( 4 . 26 ) 
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Figure 4.2. Eigenvectors for a general time dispersive uniaxial 
medium. 
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H 3 H+=rf l ' = 


C 2 

€v 

CC 


f 7 ? cc 

■rf Tft, 

vt ? 


(4.27) 


2 2 

By considering the known inverse 1 transforms of p and k and the proper 


change of variables, one can determine the function G^, i.e. , the inverse 
transform oi 


>f (n't N H ) 1 . The functions, G., that result are, 
i i i 


G < ir)= "27rN- l,ogP)8(z) 


(4.28) 


G 2 (t) = N q N| G ( #G 3 


(4.29) 


G 3U ) ■ T, T2 — 

47r(TP+Cz 2 ) l/2 


(4.30) 


where 


P = (x 2 +y 2 )' /2 


(4.31) 


2 2-1 

Similarly, by considering the known inverse transform of (k - k^) and the 


proper change of variables, one can determine the function Q^, i.e., t 


he 


-1 


inverse transform of S The functions Q , so determined are. 

i l ’ ? 


Q,= 


'/? 3 /2 2 2 . *2 

k 0 k, V£ iyw * k iV 1 

47T (K 0 N, P Z * K, N,z 2 ) ! ' 2 


477 " 


(4.32) 
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Q 2 


3/2 V2 

K, N 0 N, 

47T 


2 2 ^2 
e -jko(K,N 0 P^K,N,z ) 

(K,N 0 P 2+ K | N,zV /2 


Q3 = - k 0 8 ( r ) 

With the above functions it is possible to give an expression for 
Green’s matrix, r, and the magnetic field. 

3 

r (7,0)) = -jOU € 0 X Qj *Gj *H, (V)H*(V)8(7) 

“ i = l 


h (7, oj) = r (7, cu ) * M m (7 ) 


4.1 Field of an Electric Dipole with a Longitudinal Orientation 

The Fourier transform of source M (r) for the electric dipole 

m 

A -*> 1 A 

J (r) = 6(r) z. is M (k) = -7- r- k x z. The source M (k) is ortho- 

e 7 m (K cj £ ) m 

o o 

gonal to both eigenvectors and H^. Therefore, 


and 


B( r ) = -i(jJfL 0 € 0 0, ( r) * M m ( r ) 


B(7) = -jO/ Mo € 0 Q, (7 * * ( " j ojg^K^ ) V X 8 (7 ) Z 


(4.33) 

(4.34) 
the 

(4.35) 

(4.36) 


(4.37) 

(4.38) 


(4.39) 


79 


B(7) = l 0 ^ 0 Vxq,(T) z 


(4.40) 


or 


b(7) 



2 2 l/o 

-jko(K 0 N,P + K | N , Z ) Z 

V 

(K 0 N , P 2 + K, N, Z 2 ) l/2 


A 

z 


(4.41) 


Because B(r) is transverse to the magnetic field and the particular form of 
the permeability matrix, the magnetic field intensity is. 


h( 7) = -d— Vxq,( 7) z 

*<>N| 


(4.42) 


It then follows that. 


d(7,oj)= -1- Vxh(7)- 8(7)z 

)(x) 


(4.43) 


Notice that only one characteristic field, and accordingly one sheet of the 
dispersion surface, is involved in the source field. 

4 .2 Field of an Electric Dipole With a Transverse Orientation 

The electric dipole is directed along the y axis in order to be able to 
check Clemmow's results. The source is M^(r) = - ( J V x K 6 (r) y. 

-• / i \ °-» A 

In the transform domain this becomes M (k) = + 


(= 4 ) 


k x y. 


M m ,( k) • NT V* k x £ ( k x S ) + ( ) k x y 


(4.44) 
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( 4 . 57 ) 
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GO 


Q, * G, = (27T) 


-3 


III 


e Jk ’ r d £ d Tj d£ 


= (27T) 


-3 


t" N l ^ 2 (L 0 T l f 2 *L,T l C 2 -ko) 

JJJ e'^ “* PdPd^'d^ 


= (27 r) 


O 0 -CD 


CD CD 

*<If 


N|7> 2 [l 0 T,/> 2 + L,T, (f 
J 0 {PP)e~ l ^ z dPd£ 


-4 


= (27rr 2 (K,T,) 


-• ^[L 0 T,/» 2 + L,T,£ 2 -k 2 ] 

'* J J dPdi' 


O -CD 


„ (K|T,) I/2 f J 0 KPP)e 


47T 


I 


( , [c' 2 -L 0 T l /> 2 -k|] 
-a/l 0 t,/> 2 -ic§ I a/kT^i * I 


d/> 


/L 0 T,A> 2 -k 2 


dQ,*G, -(K,T,) ,/2 y 

ay = 


CD 


47T 


if 

0 


J,l/>P)e 


- L 0 T \P 2 “ *0 V^ | N I 2 


d/> 


L 0 T,/> 2 -kn 


(3 Q| * G j K j y 


GO 


21 


Integrate by parts 


J J,(PP)e 

o 


-\/l 0 V> 2 -kg |v^N| z 


dp 


u = 2 1 

-Iv^s; 7 1 e- 1 


L 0 T,P 2 -ko 


(4.58) 


(4.59) 


(4.60) 


(4.61) 


(4.62) 


(4.63) 


(4.64) 



V = J, (PPidP 


* 2 Qi»G, 

dz <3y 


WT' S 9". I) 



-17m 


|2 


® - /l T - 1{2 

LqT, J* J^(/>P)e V 0 1 0 


PdP 


} 


L 0 T, P ~ k 0 


(4.65 


K, y f e z I -|yK,N l Z | 

47f T (S 9"- Z) 1 P — 1 



The integral of Equation (4.66) is Summerfeld's formula (page 34, Functions 

12 

of Mathematical Physics , Magnus and Oberheit inger ), 




if 


r 



4.6 
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d Q|*Gj 
dz dy 


q 2 * g, = 


e iktT d £ d Tf d£ 


-« N,/> (L,T 0 /> + L | T 


.£ 2 -ko) 


= (27D 


? y? .-‘fr 

0 0 -CD N, /> 2 L|To/> 2 +L,T, C 2 -ko 


= (2 77") 


-2 . - 


*7/ 


JqC^P) e J?,z dPd£ 


0-® p L ,T, /> 2 + L, T| ^ 2 -ko 


= (27T)” 2 CK,T,} /2 


rr J 0 (/>P)e i * v ' K| ' N| 2 df d£' 

J J / 2 • 2 2 

0 -® + L ,T 0 P -k 0 


(K,T,) I/2 f J o (/ , P)e' v/L,T ° /,2 " k 0 i VkTn, z I d 


J* _£b( 


L I T 0^ * k 0 


dQ 2 *G, 


(K,T,) 


J,(/>P)e 


T -It 2 

L ro r *o 


K i N i z l d p 


r 2 2 

LiT^P ~k n 


do 2 *G, k, x , f -y Ll T„/>*-i,* 

a 7a » " ~ 4 tt -p (5 a n 21 J v 10 0 
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Integrate by parts. 



z 


du 





LiT 0 PdP 


L,T 


-k 


2 

0 


dv = J,(/>P)d/> 


-J 0 {PP) 




(4.77) 


The integral of Equation (4.77) is Sommerfeld's formula (page 34, Functions of 

12 

Mathematical Physics , Magnus and Oberheitinger ). 


d 2 Q s *G, k, * 

did, " 47T "p" 


(sgn. z) i 


-jkolv^i 


-|Vk^z| e' ik °V^ 


jko A /K iN 0 p 2 + k,n,z 2 


K,N 0 P 2 + K( N, z 2 


1 

I 


(yK,N 0 p and K , N ( z rea!, < arg ° / - k o ^ ] 


- a f <3 2 Q 2 *G 


B2xU)= K, dx \ dzd* 


Bg y ( r ) = 


* oaf ^ 2 Q2 ^g ( 


k i (3 y \ <3 z < 3 x 


B Zz( r) 'K|N 0 / x Q 2 (r) 


B,(T)^(s 9 n. z.VxXje-^^'l-lv^; *|^} 


(4.78) 


(4.79) 


( 4 . 80 ) 


( 4 . 81 ) 


(4.82) 


where 
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7 ) = 

47 T 




Iv^Mh z l(^) ^ 


y \ / j^o K oN| y K 0 N, y 


jk 0 R i 


M 0 y z f-(x 2 -y 2 ) e ' jk ° R| jko K o N iy 2 / j \ e ' ik ° 

Butn^AS; pHr--HT“— Rl — O-i^R,)— rT 


p“ L p 




+ 4 ¥ (Sgn - Z) ~ 


(x -y ) - jk0 |v^ 


o #-*\ _/X ° y / x <3 / I \ f - jK 0 | K . N, 

B, y (r)= W (sgn. z)-(-)(e 


-JkoRi 


K,N, z — ~ 


i^L Ugn . 2) /2«. ^-jkoIA^ z Li./^Tn^ ,\£ 


m 


jk 0 R| 


K. N, z 


| Vk|N| z | / jk 0 K 0 N| x K 0 N| x \ -jk 0 R| 


R I R? / 6 


Moxyz v/KiNi rr 2 jk 0 K 0 N, / j x -| e " ik O R i 


{[* 


47T p 2 L L p2 Ri 


WV R. 
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2Mo*yUgn. z) -jk 0 1 VkT^ * 


(4.89) 


47TP 


B, 2 (r) = 0 




-jk 0 R z 'i 
R 2 J 


(4.90) 


(4.91) 


where 


2 2 

R z = (K|Nq P *KjN|Z ) 


(4.92) 


B 2 ,(T)=^(sgn. 2 ){^[e- ik ol^ z| -|v^ 


-]k 0 R z -I 

Z 1 R 2 J 


I z | -v (M 2i . w\- w , j 

P \ R 2 R 2 / J 


(4.93) 


^0 


47 r " v K,N| 


z f lx 2 - y 2 ) e~ ik o R 2 jkpK^ox 2 / j \ e "i k 0 R 2 ^ 

p2 \ p2 R 2 Ro V kn^2 R? J 


. ll ,x2 - y2) -|kolv / M^ 

47f tS9n - 2)— ^r- e 


(4.94) 




■jko R 2^ 


(4.95) 
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Because of the form of this particular source, that is, a source that is 
orthogonal to the third characteristic magnetic field, H 3 (k), there will be no 
component of the field parallel to the third characteristic field, B 3 (k). 

This will be true of any electric source since then the Fourier transform of 
the equivalent source M m <k,co) will be transverse to the wave vector k. 

To arrive at the total magnetic field, B(r), the components parallel to 


each characteristic field must be added. 
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^ 3 

BU) = X B 


i = I 



( 4 . 100 ) 


-* M 0 V K | N | z 


B- ( r ) = 


477 - p 2 


-(x 2 -y 2 ) ( e ~J k o R i e - i k o R 2 


a -j R o R 


ik 0 K 0 N ,» 2 (l-j^ ) ^- 1 - *jk 0 K,N 0 x 2 


i X ,-»‘o R «! 


( 4 . 101 ) 


r 2 


- MoV^jN, xyz 


By (r ) = 


47 TP 


-2 I e‘ jk o R i e' ik o R 2 


/ j \ e' ik ° R| / j \ e ' ik ° R2 

- ik o K » N ' ( "vk ) — *i k oK,N 0 ( 1-rk) 


( 4 . 102 ) 


' R 7 


^0^2 ' Rg 


~ j k o/*o k i 2 * 


B,(r ) = 


(l-— ^ 
V k ft R , J 


47 T 


O n 2 


-jk© R 2 

r \ 


( 4 . 103 ) 


To the author's knowledge these results for a general uniaxial medium 
of the type considered here have never been derived before. As they stand, 
they represent the Fourier time transform of a dipole field with a time de- 
pendence of a Dirac delta function in a time-dispersive medium; equivalently, 
they represent the space-time solution to a time harmonic source. Several 
remarks should be made concerning the above example. First, each component 
parallel to the characteristic fields is significantly different from the 
total field. Note that the terms with the exp z !) dependence are 
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not even evident in the total field. Although the interference phenomenon 
is well known, it is particularly emphasized in the above. Observe terms 
like 

/ 'i k O R ! -) k o R 2 \ 

[*-*r-*-*r) 

As R approaches R the beating phenomena becomes less rapid. R may 
approach R in either one of two ways. R may approach R along certain 
directions in space or the medium may be degenerate. When the medium is 
degenerate two of the eigenvalues are equal and the corresponding sheets of 
the dispersion surface are the same. For this uniaxial case the medium is 
degenerate when K /K equals N 1 /N q . Thus, it is immediately apparent that 
free space is a degenerate medium and for this reason there is a certain 
amount of arbitrariness in the characteristic fields. 

For the moment, let N = N 2 = 1 in order that a comparison can be made 
with the work of Clemmow. P. C. Clemmow published a paper entitled, "The 
Theory of Electromagnetic Waves in a Simple Anisotropic Medium," in the 
Proceedings of the IEE , Vol. 110, No, 1, January 1963. In this paper, he 
gives a method to find the exact fields due to a time-harmonic source in a 
nonspace-dispersive uniaxial medium. The above derived expressions for 
N = I compare exactly with those of Clemmow. He shows that such a field is 
related by a scaling procedure to a corresponding vacuum field. The vacuum 
field is expressed as a superposition of a transverse magnetic field, J6 
which the magnetic vector is everywhere perpendicular to the axis of symmetry 
of the anisotropic medium and a coplanar transverse electric field; and dif- 
ferent scaling is applied separately to each partial field. Only in passing 
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is it mentioned that these scaled transverse fields are two of the three 
characteristic fields in the uniaxial medium. Thus, the full signficance of 
the characteristic fields is not utilized. Since only two characteristic 
fields which cannot span the three-space are used, it seems questionable on 
the outset whether the field due to an arbitrary source (both electric and 
magnetic) can be described only in terms of these. However, the difficulty, 
which Clemmow does not mention, can be indirectly circumvented by using a 
superposition of the fields for each type of source. Another difficulty, 
the task of resolving the vacuum fields due to an arbitrary source into TM 
and IE fields, is not a trivial matter. The whole procedure depends upon 
the ability to split the source into components that excite each type of 
field. 

In l ight of the proceeding discussion let us go back to Equal ion (1 . b) and 
observe tie explicit role played by the characteristic fields. We wish to 
show that by the use of the characteristic fields the source resolution is 
automatic and is not simply an artiface for a particular problem. Let us 
begin by saying that the source field H(k) may be expressed as the weighted 
sum of the characteristic fields. Thus, in general, all three characteristic 
fields will be involved, although for some source distributions the weighting 
factor may be zero. Equat ion (4.4) now becomes, 

3 

Z GiGHHj* — )(x) € 0 M m (4.104) 

i= t 

But G^jHj = SjNHj =/* 0 SjBj . Hence, we have 

3 

Z ^ « S j B i = j € 0 M m ( k ) 

i = I 


(4. 105) 
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The weighting factors a. are found by multiplying Equation (4.105) on the 

J 

left by Ht\ 

o 



(4. 106) 


At this point the prominaht role played by orthogonality is easily seen. Also, 

the characteristic field will not contribute to the total field if and only 

if hT M (IT) equals zero. This is exactly the case for H and an electric 
J m 3 

source current, since HL equals k and M (S’) is transverse to k. Since 

o m 

N = I, H = kx£ and H = k x k x £, then H (r) and H (r) correspond to the 
scaled TM and TE fields, respectively. The important factor to be remembered 
is that the TM and TE decomposition applies only because the two character- 
istic fields are TM and TE. However, the characteristic field decomposition 
(spectral decomposition) applies for any medium. 

Let us go one step farther to produce another result. As was previously 
stated, Equations (4.101), (4.102) and (4.103) may be interpret ated as the 
time Fourier transform field of a source J( T, t) = 8(x ,y, 2, t) y. Now assume 
the constitutive relationship is independent of k and go. This assumption 
makes the Fourier inversion possible. The results are almost immediate. 


B x (r,t) = 


Vkx 1 


47 r 


{ 


- ( x 2 — y 2 ) / 8^ c c 

D 2 \ R rT 


(4. 107) 


+ K 0 N, 


r: 




-r R i>+ 


S(t~F R |)\ K,N 




S(t c ^2 


)} 
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B, 


r „ Hi f -2 ( SI'-f R,) 

' ’ 47 7* P 2 l"p 2 ' R; R 2 / 


(4.108) 




B. 


. ~M° k/ 2 No^ (1 S '(,-J- R ) + S«-c r ;> x 

47TR? xc c 2 R, • 


(4.109) 


The support of the fields is two ellipsoids expanding with time in contrast to 
the support in free space which is one sphere expanding with time. Thus, an 
observer at an arbitrary point in space will be cognizant of two wave fronts 
not one. Although these ellipsoidal wave fronts are not of the same shape as 
the sheets of the index dispersion surface, they are related. The ellipsoidal 
surfaces are of course expanded or contracted ray surfaces, which may be given 
in parametric form by 




VkS 

T-V k s 


(4.110) 


A problate ellipsoidal index surface will have oblate ellipsoidal wave fronts, 
and vice versa. Another point is made apparent by Equations (4.107), (4.108) 
and (4.109). No point on the wave front can travel faster than the velocity 
of light in a vacuum. Therefore, restrictions are placed on the arguments 
of the Dirac delta functions which in turn places restrictions on the 



permeability and permittivity matrices. Not only must the permeability and 
permittivity matrices be positive definite but also 1/N q , l/N^ < and 
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5. SPECTRUM OF CHARACTERISTIC WAVES IN AN ISOTROPIC COMPRESSIBLE PLASMA 

Very few problems are solvable in closed form. As has been shown, the 
dipole in a nonspace-dispersive uniaxial medium is solvable. No doubt even 
some relaxation of the nonspace-dispersive restriction can be made. The 
biaxial and magneto-ionic problems, however, as yet have not been solved. 
Numerous attempts have been made, as the literature will testify, but all 
results have involved at least one unevaluated integral. The closed solu- 
tions for the latter two media have not been obtained even though they are 
nonspace-dispersive. The problem of a dipole in an isotropic compressible 
plasma of N mobile ion species is solvable in closed form even though the 
permittivity matrix is both nondiagonal and is space-dispersive. The reason 
why some problems are solvable and others are not seems to lie in the form 
of the determinant al equation or equivalently the characteristic equation. 

The eigenvalues of the solvable problems do not involve radicals whereas the 
eigenvalues of the unsolvable problems do. These observations should be more 
apparent after a comparison of several problems is made. For this reason and 
to further illustrate the usefulness of the spectral decomposition, it is 
instructive to obtain the solution of a dipole in an isotropic compressible 
plasma of N mobile ion species by the three-vector method. 

First we must show that the permittivity matrix for the lossless plasma 

is Hermitian. The inclusion of a static magnetic field B will be made since 

o 

the degree of difficulity for proof is not increased. 

Lemma 5.1: The permittivity matrix of a lossless compressible plasma with a 


static magnetic field is Hermitian. 
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Proof: Assuming the fluid model for a plasma is applicable the force equation 

and the continuity equation can be written for each species of mobile i<i>n 
(including electron). 

m j -|j L - q j$+ ( m i a -/N J ) V n j-»- q j B XW'j = 0 t 5 - 1 * 


cm-, _ 

-yp-+NjV^ = 0 (5.2) 

t h 

where a is the speed of sound for the j ion given by 

j 

a 2 = Y.K T ,/m. 

J J J J 

v = ratio of specific heats at constant pressure to that at constant 

j 

volume 


K = Boltzmanns constant 
Tj = ion temperature 
m = ion mass 

S 


Equations (5.1) and (5.2) with Ampere’s equation is sufficient to derive the 
permittivity matrix. Eliminating the density n^ from Equations (5.1) and 
(5.2) yields 


aVj 

W 


of VV-^j + ^j 


M. 

3t mi at 


(5.3) 


where 




B 0 X 


(5.4) 


The temporal and spatial Fourier transforms of Ampere’s law and Equation (5.3) 


ire, respectively. 


-jkxH(k,OJ) = jOJ€ 0 E(k,GL>) + X Qj n j v j + 

q* _ 

( Cx ) 2 - O f k k T - jOJil j ) Vj ( k , CJ) = - jcu — E ( k , U) 

- Qi - 

R j Vj (k,GU) = - jOJ -j^r E(k,CJ) 

Therefore , 

Vj(k,(J) = -jCU-^rRj'EU,^) 

2 

- jkxH lk,GJ) = jCL >€ 0 Jl” X m .^ o J + 


From Equation (5.9) we see that the permittivity matrix is 




where 
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The index j is summed over all species of compressible ions. As seen from 
Equations (5.7) and (5.8) and the fact that is skew-Hermitian for real k 
and «, Rj is Hermitian. Hence, the inverse of R . is Hermitian. It then fol- 
lows that the permittivity is a Hermitian matrix, i.e., 

(5.12) 

This derivation has assumed that collisions are negligible (lossless) and that 
the ion pressure for each species obey separate adiabatic relations. 

From the previous derivation for the permittivity matrix it is observed 
that, 


2r*T, 


(OJ - Qj k k ) 


(5.13) 


and then 



I - 


2 T r T 

o j k k 

" 2,2 T. 

(a j k -oj ) 


Therefore, the normalized permittivity matrix is 


(5.14) 


K = 



I ♦ 


I 


2 ,2 

a ■ OJ ■ 

J J 


2 2 2 2 
j CU (o j k -OJ ) 



(5.15) 


where u. is the angular plasma frequency of the j th species, expressed by 


I rrijf: 


-5.1 Eigenvalues for an Isotropic Compressible Plasma 

Since Gg(k,w, \) and Gjj(k, co, \) have the same eigenvalues, it is sufficient 

and easier to evaluate them explicitly through det G (k,w \) =0. It will 

H 
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become apparent that the form of K is the simplifying factor, £ is of the 


T -1 -1 

form bk k + a, and K is of the form K = 1/a 

-i *r 


I - (k + a/b) 


-1 


1 -*-*T 
k k 


G (k «,X) = (-k x K kx » X) . Therefore, the second term of K has no 
H ~ — 


effect upon G , i.e., 
H 


G h Co;,X ) = (- -i— ITx kx - x) = - -i- [tr»T— Ik-oX)] 


(5.16) 


Then 


k 2 2 

det G H (k t u»,X) = — X(X-~-) =0 


GO 


(5.17) 


From the previous section it is found that a = I 1 ~ £ Hjj- I and 

\ J GO 

2 2 
a . go . 

b = [ £ - - — | * The eigenvalues are independent of the com- 

2 y 2 _2 2 . 
go (a k - go ) 

j _ 

plicated term b. Since a is independent of k, both eigenvalues X^ and X^ 


(f <1 


Sl = S 2 = (k 2 /a - k 2 ) 


which are degenerate will have only one sheet^ 

1/2 

sheet will be propagating for w > f 2- wj and nonpropagating for 


The 


1/2 


“"(H) 

In summary, there is only one degenerate transverse sheet due to the 
eigenvalues for an istropic compressible plasma. 

5.2 Transverse Part of the Field Due to an Electric Dipole 

Since the eigenvalues X and X are degenerate, arbitrarily choose, 

1 2 


E,= kxz , E 2 =kxkxz, E 3 =k 


(5.18) 
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Also since the plasma is isotropic, arbitrarily choose the electric dipole 
source orientated in the z direction. A magnetic source cannot excite plasma 
waves; therefore, at this time we are not as interested in the fields due to 
it. With this choice of eigenvectors, does not enter into the picture. 


E+ K E 2 = (!TxTxz) + (aI + bR J J )(k*kSz)-- 0 p 2 k‘ 


(5.19) 


E 2 (iT,aj) = k 2 nTxlTx^) 


(5.20) 


By partial fractions. 




(5.21) 


Therefore. 


-jVok.r 




(5,22) 


5 .3 Plasma Waves 


The eigenvalue for the longitudinal or plasma waves is zero. This implies 
that propagating plasma wave must come from the real zeros of K . Since 

g 3 mT ’ 






(5.23) 
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After obtaining a common denominator and simplification, one obtains 



x 2 t 

Except for the factor k in Equation (5.24) the numerator of E^ K E g is a 

2 j . 2 

polynomial in k of order N. Therefore, E^ K E (discounting the k ) will have 

2 

N zeros. Since the numerator of Equation (5.24) is a polynomial in k with 

2 

real coefficients, the roots of k will occur in conjugate pairs. Hence, it 
is possible if all of the roots are real to have N sheets to contribute to 
plasma waves. Restated, the maximum number of real sheets (that cause propa- 
gating waves) for an isotropic compressible plasma is one sheet per compres- 
sible ion species. This is further verified by obtaining the determinant of 
K and comparing the zeros with the zeros of E^ K Eg. 


del 




2 , . 2,2 

OjOJjk 


l—j 2 2 2 2 

j oHojk -of), 



(5.25) 


The zeros of E* K E are the zeros of 
3 3 



Equation (5.26) can be put into a form 


(5.26) 


7T ( k 2 — k 2 t ) = 0 

} 


(5.27) 
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where k is the zeros of (5.26). The product of the zeros is 
J 


2N 

_ L 2 

7T kj - 2 

J 7T a-, 


-e4 

j cu 


The sum of the zeros is 




(cu 2 -w- ) 


2 2 2 

From Equation (5.28) one observes that for co < 2 to , IT k < 0, implying 

j J j J 

that there will be at least one negative zero, i.e., at least one nonpropa- 

2 9 

gating plasma sheet. For to > 2 to" if N (number of compressible ion 

j J 

species) is odd, there is at least one propagating sheet. 

Now consider the longitudinal fields, 


ieJke 3 )" 


_ / . 2 (Jj 

TT I k : 


. 2 ^ .2 0J 

k TT I k ? 






j a-, ' V a 


k 2 IT 


\ a j / J a i \ a * / J 


+ -1 2 

l/a is the residue of (Eg K Eg) at k =0; therefore, the second term of 

2 

Equation (5.30) does not have a pole at k =0. 


E 3 lk,CJ) = -jaJ Mo (-k‘ 0 2 )(EjKE 3 r'E3EjM e (k,a)) 
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Mgl k ,(JJ) - Z (5. 32) 

Let 

E 3 (k t CJ) = E 3 , (k,CU) + E3 / (k ,0J) (5.33) 


where 


e 3 u,o)) 



(5.34) 



Since the medium is isotropic, the integral is independent of 9 and <j>. 
Therefore, for the sake of computation, let 0 = 0 and integrate over <$>' 



The integral should be integrated over the closed contour, C, which is a 
path along the real k-axis and an arbitrariy large semicircle in the lower 
half k-plane. The path should be such that poles on the positive real axis 
are included in the contour, but poles on the negative real axis are ex- 


cluded. Denote all of the zeros of 



on 


1 to N. 


the positive real axis and in the lower half plane by k. for j ~ 

J N 2 2 

Then the denomonator of the integral of Equation (5.36) is k n (k k ) # 

j J 

Then Equation (5.36) becomes, 


VV T z N 


Ejlr.tt))* 2 
a k o 


4777 


I 

u-l 


. 2 O) \ V W i (j)/, 2 Of 
(0-1)77- k u --2 )-L-iT (k u — g 


J 


'J 


i a i ' 


a; 


k u 7T lu) (k u -k j )7T(k u + kj ) 


(5.: 


ri k u r 


if all of the poles are simple. In general, 


£,'(7.0)) = 


ia, t^yv T -l £ 

ak 2 27rr c 


Res 


2 « 

. 2 CU 2 \ U) / 2 C*J 

(a- 1 ) 7T (k - — (-Z.— V ( k - — 


J 




j aj 




ri kr 


(5. 


5.4 Total Field Due to an Electric Dipole 


e(7,g o) = e 2 (7,oj) + E 3 (7,o)) ♦ E 3 (7,oi) 




The residues inside C are on the positive real axis and in the lower half 
k-plane . 

Note that for N = 1 this result checks with Equation (28) of Hessel and 
13 

Shmoys, "Excitation of Plasma Waves by a Dipole in a Homogeneous Isotropic 
Plasma,” Proceedings of the Symposium on Electromagnetics and Fluid Dynamics 
of Gaseous Plasma , Microwave Research Institute Symposia Series, Volume XI, 
Polytechnic Press, 1962, pp. 173. The "modal decomposition” of their paper 
is not the same as decomposing the field along its eigenvectors; however, the 
relationship between the two can easily be seen. The relationship is, 

EjrTcj) =E 2 (f\CJ) + E3(r|a;) (5.42) 


VjL/ n 


(5.43) 


A ^ yy 

3S(r,u>) is not entirely longitudinal or transverse, but <7 P 00 E (i^w) is 

3 


longitudinal . 
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5 . 5 Pressure in an Isotropic Compressible Plasma 


th 


The partial pressure variation for the j component is 


Pj = a* m j"j = y j* T i n J 


However, from Equations (5.2) and (5.8) one finds that 


n,(k,a» • r-'eu.oj) 


Using Rj 1 given in Equation (5.14), the partial pressure P^(k,u0 

_ JO; q i N 1 _ T _ 

P,(k,(D) = ' ' T k 

(aj k -Ct) ) 

whereas the inverse transform of P^.(k,o.) is 

Pjlr.GJ) = “N j q ] (- ! 4^-)* vTE(r * aJ) 


Thus, the total pressure variation P(r,w) is 


(Jj 

e"V 


7 t _ * 


PU.GU) = ZPj^,OJ) = -\Z N iQj— 5: ^Fr)* V E(r » w) 


Notice that only the longitudinal field contributes to the P r(: 


(5.44) 


(5.45) 


( 5 . 46 ) 


( 5 . 4 ? ) 


(5.48) 


107 


6. RADIATION FIELD OF AN ARBITRARY SOURCE IN 
A LOSSLESS LINEAR PASSIVE MEDIUM 


6.1 Introduction 

For the case of an isotropic medium, the radiation field for an arbitrary 

current source is well known. The problem of finding the radiation field of an 

antenna in a cold magnetoplasma such as an ionized gas in a constant magnetic 

14 15-19 

field has been solved by Bunkin and a number of workers. However, such 

solutions have been limited to a particular nonspace-dispersive media. It is 
desirable to determine the radiation field for a general lossless linear space 
and time-dispersive medium. It is preferred to find the general field solution 
for any zone; however, as yet no one has achieved this for the cold magneto- 
plasma or even the biaxial medium. For a cold magnetoplasma with one mobile 
charged particle species, the problem is to solve nine first-order partial 
differential equations. If a harmonic field is assumed, six of the unknown 
variables may be eliminated to yield three second-order partial differential 
equations in say the field variable . Usually the system of three equations 
is attempted to be solved by the method of Fourier transforms. The solution 
can then be expressed as a volume integral in Fourier space. Thus, the dif- 
ficulty is a triple integral of a function with a complicated singularity. 
Invariably, attempts at such a solution are expressed in at least one unevalua- 
ted integral. 

With the results expressed in terms of an integral , it ..is difficult to 
make comparisons between fields of different sources and to interpret the 
physical processes. Thus, it is desirable, and for some purposes sufficient, 
to find an asymptotic solution for the radiation zone. Notably, two asymptotic 
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integral methods have been used in the past, the steepest descent or saddle 
point method and the stationary phase method. To first-order, they produce 
essentially the same results. For our purposes the stationary phase method 
will be used since it will be seen to yield many physically interpretable 
results. Only lossless media will be considered throughout. This is not; a 
severe restriction since in a lossy medium the concept of a radiation field 
is not very significant. 

First, the principle of stationary phase is applied to a very general 
linear system which encompasses both space and time dispersive media. Then 
the results are particularized to some systems including warm plasmas, which 
are derivable from dynamical models. Notable physical interpretation of the 
mathematical results are made. 

6.2 Stationary Phase Method for Arbitrary N-Vector System 

Consider a system of equations whose Fourier transform is 

0 - j OJ U 

or 

F - C (6.2) 

where the N-vectors C and F are the source and field vectors, respectively, 
t h 

Also, the IT order square matrices 0 and U are skew-Hermit ian and Hermitian, 
respectively. All of the quantities may be a function of all the transform 
variables. Then similarly to Section 3 an eigenvalue equation may be defined. 


F= C 


( 6 . 1 ) 


o Fj = v \ y Fj 


, i = 1 , 2 , • • * , N 


(6.3) 
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Analogous results concerning the eigenvalues and eigenvectors may be inferred. 
Thus , 

N N 

F(k ,(j) = X Fi (k ,oj) = y“'Z s "i' r i c (6 * 4) 

i= I i=l 

where I f = (F* U F. ) _1 U F. fT and S. = v. - jco. In a form that we will use, 
1 i = 1 = ii l l 

the portion of the field that is parallel to the characteristic field F^, is 


F j (k,a>) - (Ff^Fj f 1 Fj f|c 


( 6 . 5 ) 


with the normalization of F^ such that the components of F^ contain no singu- 
larities. Therefore, the inverse Fourier transform is 


y tft) = (27 t) 


-4 


fffj 


- 00 


(F; + ^F; f'F, 


F: + C 


^('oJt-k-r ) d 3 kdaj 


( 6 . 6 ) 


v 


The field solution as described above is not unique in that an arbitrary source 
free solution may be added to it. Furthermore, we wish the solution to describe 
the physical model. Thus, all waves must originate at the source and no source 
free solutions are permitted. In order to obtain such a unique field a "radia- 
tion condition" must be enforced upon the class of possible solutions. Now the 
solution is a frequency spectrum of time harmonic waves as the outer integral 
of Equation (6.6) stipulates. Let it be required that each of these component 
harmonic waves obey the "radiation condition. And let the Fourier frequency 
variable, o>, be slightly complex, i.e. , replace jco by s = a + j<*> where a is 
small and positive. Denote the frequency component field by t, . 
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GO 

f; (s,T) = (277-) 3 jfjf^(F^ Tf\ F 1 T , F i F'fce (s, ' ik ' r ^ ) d 3 k (6.7) 

-00 

With this, after the evaluation of f^s^r), the required field for the physical 
problem is, 

00 

2(rl\)=(2TrV fU^ f i (s ’ r )] dcu (6 8) 

J Ls-moj J 

-00 J 


The reasoning for the mathematical manipulation of the Fourier frequency variable 
is as follows. With O positive the harmonic frequency component f^(s >■*> is 
increasing exponentially with time. Now if a source-free wave of finite ampli- 
tude is propagation from infinity then by the time it reaches a finite distance 
from the source it will be small compared to the exponentially increasing com- 
ponent f^Cs^r). Thus, if only the exponentially increasing waves are sought, 
the source-free waves will be omitted. Taking the limit finally produces the 
Fourier frequency component satisfying the "radiation condition." Let us now 
asymptotically evaluate f^Sjr) for large r. 

The denominator (F^~ 7^ F^) is changed with the introduction of the new 

variable s, and its new form is approximately, 


Ft ^ F; (s,7)» F+ 77}Fi {So J o ) +ds A 



+d rV y 




( 6 . 9 ) 


Ill 


If (S o ,Y o ) is a root > i.e., 


’Ity* i] <S o. 7 „)=°- 


then for 




(s,Y) 


to be zero when s = s q + ds requires that 


y = y 0 -d/ = y 0 -ds 


A 

,ds 


[f^Fi] (s,?)r/r V^F^Jts.r) 


^ ( 6 . 10 ) 
sqi r 0 


Let s - ju), Y = jk and ds = (T„ Then the new root is shifted to 


= a + jo), r=\k-a (-&- [pf^ Fj] (OJ,k) r/r V- [r^Fj] (0J,k) 


(6.11) 


Since the operator Tfj is skew-Hermit i an, the new root Y has been shifted off 
the imaginary axis. With this in mind, perform the first integration of f (s,r) 

♦ A 

on a variable k| ( = k«r which is parallel to r with the transverse variables kj^ 

held fixed. Evaluate the integral by the contour integration method. By the 
Cauchy residue theorem we have, 

jR 


f {Ft 77/ Fj f'Fj Fj + Ce (st ‘ r " r) dy„ 
-j R 


+ J (Fj + T7J Fj ) 1 Fj Fj + Ce (st y,|r) dr n 

L 


( 6 . 12 ) 


= -27Tj^ residues of poles in the contour 


where L is a semicircular path of radius R in the right half plane in the clock- 
wise direction. Then by substituting — Y || for p and r for t in Lemma II of 

20 

Transform Calculus by E, J. Scotty we have 


Mrrv /(F i t ^F i )' l F i Fi t Ce ( * , " > " r, dy, l = 0 (r>0) 

R-*ao l 

provided the components of | (F* 7fJ F^ -1 F F* c e St | are less than M/ 
when | Y|| | > R where M, K are constants and K > 0. Therefore, 


( 6 . 13 ) 



K 
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J&??u I 

R -+■ 00 * 

-jR 


CO 


'F i r l F i F^Ce ls, - > '" rl dr ll 


277-iZ residues in right half X M plane 


.-I _+ _ (st - jknr ) 

) Fj Fj Ce 11 dk. 


= -277'X residues in right half /„ plane 


(6.14) 


(6. 15) 


But it was previously noted that the poles in the right half plane are the poles 


such that r 


■v- Fj^Fj (O),io)<0 


>for <7 > 0. Since the source, 


C, is assumed to contain no singularities, then for simple poles the residue at 


a pole is 


-I _ r-t ^ _ (st - X- r ) 


(r-V^F^Fj f Fj Fj f Ce 


Thus , we have 


(6.16) 


^Wfj(s,T) = i( 2 V)~ Z JJ fr-V k Fj^Fj Fj F i 1 Ce j(<J,_k ' r) d 2 k 


(6.17) 


ire 2 i+ is the portion of the dispersion surface in which 17^ F . (w,k) = 0 

' fr-V- ff^F, (OJ.k)/-^ [Ft^Fj] (CU.tT)) < 0 


Now perform the last two integrals by the method of stationary phase. Hence, 

the major contribution to the integral for large r is from those points on 

-► -+■ 

S. where the exponent k.r is stationary. These are points where V T k* r = 0 
i k 

or, equivalently, points where the normal to 2 i+ is parallel to r. In the 

vininity of a stationary point k we can express k.r as a second-order surface 

K 

in the variables transverse to r. Let the transverse variables k , and k , „ 

_L1 -12 
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be in the principal directions. Because k.r is stationary and k ^ and k ^ 
are the principal directions, k ^ , k ^ and k ^ form an orthogonal coordinate 

system. Also, let p^ K and p be the curvatures of the k surface at the sta- 
tionary point k^ associated with the principal directions 1c ^ and k^, 
respectively. A positive value of curvature implies a concave curvature to ?, 
whereas negative implies a convex curvature to r. Then we have approximately. 


k* r = k 


K*' + [ 


2 ^ IK^ll k llK^ + 2 ^2K^ k l2“ k l2K^ r 


(6.18) 


By the principle of stationary phase the asymptotic form for r approaching 
infinity is, 

^fj(s,T)^j(27rf 2 I^V; Fj*/^ 7 Fj j (CJ,ky ' F. F^C ^t-k^r) 


a) r 

•//e-'zt 


^W k ll k llK ^ + ^2K^ k l2 


k l2K ^ ] r 


(6.19) 


dk 1) dk i 2 

where jp represents the sum of all stationary points on The integral in 

Equation (6.19) can be evaluated after possibly slight modification from most 
integral tables. Hence, 


fi(s,0«j(27rf'X (r-VHFj^Fil ) F,F*C- 

s-j U) K \ L J / r 


J " j -J(sqnP,+ sgn ) 


•A W2\ 


( 6 . 20 ) 


j(OUt-kT) 


But the product of the principal curvatures is the Gaussian or total curva- 
ture, i.e. , p p = . Therefore, the asymptotic form of the field (r,t) 
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for r approaching infinity is 


CD 


7, <r,t) * j 1 2 7 T)'*fz (?• V^F^F, 


Fj F,'c 


-00 


Vra 


• j-^ (sgn P t + sgn P z ) 


j(GJt-k*r) 
e dOJ 


(6. 21) 


And of course the total field for the system is. 


N 


t’(m) = X^t) 

i = I 


( 6 . 22 ) 


J 00 t 

For the special case of a time harmonic source ^(r,t) = (r) e ° ^ the 

Fourier transform is 


C(k,ai) =2 7TC("k)8(aj-cu 0 ) ( 6 . 23 ) 

where 6 is the Dirac delta function. With this equation. Equation (6.21) is 
easily evaluated to give 


^t)=X r-ViT F^F, 


(CU 0 , k)) FjF^CU) 




e 


-j-^(sgn P { + sgn#,) 


j((D 0 t -k-r) 


(6.24) 


6.3 Application of Section 6.2 to an Arbitrary Six-Vector Electromagnetic System 
The general description of an electrodynamic system can be expressed in 
terms of Maxwell's equations (Faraday's and Ampere's laws) and an auxiliary set 
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of equations. The auxiliary set of equations may be partial differential, 
integral as the expectation of a Boltzmann’s current, or simply an algebraic 
equation as the Lorentz equation for moving coordinates. However, the aux- 
iliary set will be coupled to Maxwell’s equations by common variables. Thus, 
often times when only the electromagnetic field is desired, it is convenient 
to eliminate all the variables except the electromagnetic field variables 


y - 




In doing so the partial differential equations usually become 
of higher order and other complications set in. Alternately, the space-time 
Fourier transform of the general linear system of equations may be taken and 
the extraneous variables eliminated by algebraic means to produce a six- vector 
equation of the form, 


OF - jCUUF+C 


(6.25) 


where 0 and F are the Fourier transforms of & = 


0 Vx 

L-Vx o 

tively. C ; is the Fourier transform of the equivalent source involving all 


and , respec- 


of the sources of the original system of equations. And the sixth-order square 
matrix U is defined to be the constitutive matrix for the relationship F^ = UF, 

P ) /¥?1 


where F^ is the Fourier transform of = 


j . Now Equation (6.25) is of 


the form Equation (6.1) of Part 6,2 if the electrodynamic system is lossless 
and the constitutive matrix, U, is Hermit i an. Thus, the results of Section 6.2 

follow. For this situation, however, as seen from Part 2.5 the quantity 

V7.* Tr- f 7n — 

J k 


F; ^ j j can be reduced further, i.e. , 


r 

VT[ F f^ij’2jP T , = 2|(P ei * I 


(6.26) 
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P . _ P . , and P are the total, electromagnetic and medium average power flux 
Ti* ei 9 mi 9 9 

vectors, respectively, for the i^* 1 characteristic field. The mathematical 

definitions of P . and P . are repeated for the sake of lucidity. 

ei mi . 

P ei = Re ( Ex H* ) (r ms) ( 6 . 27 ) 


p mi = " -y Fj T (V^ U) Fj (rms) (6.28) 


Also from Section 2.5 it was noted that P is normal to the dispersion surface 

th 

and is related to the group velocity of the i characteristic field. Hence, 
the radiation condition implies that 


A 

r 



(6.29) 


or that r*P = I P I > 0 at a stationary point. Further, normalize the 

Ti I Ti I 9 

characteristic field such that the length of its total average power flux 
vector is equal to one. Define the number h(k) to be as follows 


h(k)= 


(6.30) 


Using the above considerations the asymptotic expression for the characteristic 
source field i (r,t) for the time harmonic source as given by Equation (6.24) 
reduces to 


</; (f , t) : 


I 

K 


-j h( k K ) 


FjVlkKlFj e 


j(0) 0 t K K * r ) 


477’r 


(6.31) 
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The summation is over all stationary points on 2 + . And again the total field 
is given by 

6 

W,t). (6.32) 

i = l 

Note that this derivation is very general in that the only assumptions made 
are that the electrodynamic system is lossless and that the medium and the 
source are such that the condition for the contour integration of Part 6.2 is 
satisfied. Also, note that Equation (6.31) gives the electromagnetic field 
due to any source of the original system of equations and not just due to the 
electromagnetic sources since ^ is an equivalent source involving all of 
the sources of the system. 

— ggse of a Compres sible Plasma with N Species of Charged Particles 
Assuming the fluid model for a plasma is applicable, the linearized 
e Quations that describe the system are 




m 


(6. 33) 


VxW '- e °-£- e +Z 


j = l 


(6. 34) 


Pm j ^j = ^qj £ + P q]^ r i xB 0 ~ , j = l to N 


(6.35) 


p o]V - Sj , j= i to N 


(6. 36) 


where the quantities and their Fourier transforms are 
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(Fourier transform) 


E 

H 

J. 


m 


V j 


F. 


0 , 


5 


m 


mj 


% 

p 0) 
3 ) 


8 , 


electric field 
magnetic field 
electric source current 
magnetic source current 


. th 


mean mass density of the j charged component 
mean charge density of the j th charged com- 
ponent 

t h 

velocity field imparted to the j charged 

component by the sources 

th 

variation of the j partial pressure imparted 

by the source (normalized by the number density) 

"th 

mean partial pressure of the j component (norm. ) 

t h 

source term for the j force equation 

source term for the continuity equation 

ratio of specific heats at constant pressure 

"t h 

and constant temperature for the j charged 
component 

static magnetic field 

Equation (6.35) is the force equation and Equation (6.36) is the normalized 

continuity equation since the pressure is assumed to be proportional to the 

number density to the y th power. A number of assumptions are involved in 

J 

arriving at the above system of equations. First, the equations have been line- 
arized and hence will more accurately describe reality for small variations. 

The medium is macroscopically homogeneous and of infinite extent. A scalar 
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pressure for each charged component is assumed; hence, shear wave will not be 
evident in the solutions. Separate adiabatic conditions for each species of 
charged particles are assumed. The latter two assumptions are equivalent to 
the conditions for truncating the moments of the Boltzmann's equation. Also, 
collisions are negligible. These assumptions are not unlike those used by 
many investigators, however crude the model may be. 

Now let us show that this system can be developed into a special case of 
Part 6.3. Thus, one must find the constitutive relationship. After taking 
the Fourier transform of the system of equations, eliminate the velocity and 
pressure variation variables. In doing so, one arrives at the pair of Sector 
equations. 


j k - xH = i W kl +£ A 2 n Ri] E + [je-i£ R‘o kOn 

L n= I •* *■ n=l 

N 

-}ojy p r _i f ° 

n Tj n "J 


(6. 37) 


-j k x E = -j Uyx o H - J m 


(6. 38) 


where the matrix R is defined asR = (- go 2 p I + jwp B X+vP Te "k T ) 

n n mn qn o *n on 7 

and I is the identity matrix. From this it is easy to see that the constitu- 
tive matrix and the equivalent source are 



(6. 39) 
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C lk,tu) = 


N N 

Je ~ j ^ ^qn^n ^ - X Pqn ^n 

n = l n_l 


(6.40) 


'm 


Since R is Hermit ian for real k and u>, then U is Hermit Ian and Equations (6.31) 
n = 

and (6.32) for the radiation field apply for a time harmonic source. 

Uo for to the section on General Formulation of Spectrum of Characteristic 
Waves for a discussion of the eigenvalues and eigenvectors. 

It has been determined that the radiation field is due to the stationary 
points on the surface F^ 77^ = 0. However , this is also equal to 

Ff^Fj= SjFfyFj = U/j - jCU) fT yFj =0 (6 - 41) 


Therefore, for the longitudinal waves when vo = 0 , the quantity U f\ must 
be zero and its surface found. Another apparently added difficulty is that the 
quantity fTc^ may contain singularities whenever the matrix R^ becomes singu- 
lar. However, after a bit of algebra it can be shown that 

N N 

F?c'= Ei’j, * nfj m *£ vt F° n * X 

n=l n=1 

and thus can be made to contain no singularities with the proper normalization 
of the eigenvectors. 

Greater consideration will be made for the case N = 1 in a later section 

by the ten-vector method which is in some respects more desirable. 

Note also that when N = 1 and P = P = 0 the usual cold plasma model 

on n 

used by numerous investigators results. Observe that the constitutive matrix 


(Pni An PontOn <6 ‘ 42) 
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is not a function of the wave vector. Thus, V = P or P =0 and the normali- 

9 T e m 

zation is adjusted accordingly. This case, too, will be discussed further in a 
later section. 

6.5 Case of Compressible Plasma with N Species of Charged Particles 
(6 + 4N-Vector Method) 

Consider the same system of equations as described in Part 6.4. Naturally 
the assumptions involved in arriving at the equations are the same. Now, how- 
ever, let us operate upon the equations in a different manner to obtain the 
solution. Instead of eliminating some of the unknowns, we will keep all of the 
unknowns and seek the solution of a matrix equation of order (6 + 4N). This 
manner of seeking the solution does not seem to have been considered by inves- 
tigators before. For this reason and since there is a good correspondence 
between the mathematic and the physical behavior, the (6 + 4N)-vector method 
will be dwelt upon extensively. We wish to transform the given system of equa- 
tions into a matrix equation that will be a particular case of Part 6.2. The 
ordering of the equations and the normalization of the unknowns will be critical 
although not unique. The motivation for doing so is to satisfy the conditions 
of Part 6.2, i.e., make 0 skew-Hermit i an and U Hermitian. 

Order the equations in the following manner: 

1 . Ampere 1 s Law 

2. Faraday' s Law 

3. Force equations for particles 1 through N 

4. Continuity equations for particles 1 through N 

The continuity equations are to be in the same order as the force equations. 

In addition, place all the terms involving partial derivatives with respect to 
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time and the sources on the right side of the equal sign with the remaining 
terms on the left. Moreover, define the normalized quantities, ^UL as, 


V. n • (<?„ /r n P 0 „) ,n»i,2, • • • ,N 


(6.43) 


Then the above directions produce the 6 + 4N order matrix equation 


07-- A 

dt ~ 


( 6 . 44 ) 


which is given in Figure 6.1. The corresponding terms should be obvious. 
Notice that already we have gained a benefit, for the conservation of total 


energy is expressed as, 


where 


7'er-- y''7, * & 


-V T * T .iur T 


(6. 45) 


( 6 , 46 > 




( 6 . 47 ) 


W f k e T e • Mo t £ P mn tA>„ ♦ £ r„ p 0 „ u l 

L _.i _-i J 


( 6 . 18 ) 
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Now the Fourier transform of Equation (6.44) is 

OF “ jOJUF + C (6.50) 

where 0 is & with the del operator, V, replaced by -jit and U is 


U = 


Now it is easily verified that for a real wave vector k the matrix 0 is skew- 
Hermitian. Not only is the matrix U Hermitian and diagonal but it is also 
positive definite. Therefore, the conditions of Part 6.2 are satisfied and 
the results for the radiation field follows. Because of the symmetry of 0 and 
U, the usual type of eigenvalue problem can be defined. 

OF, - x j y F ; J = 1,2, * • • , 6 + 4N (6.52) 

Since 0 and U are not a function of co and are of order 6 + 4N, the 6 + 4N eigen- 
values will be a function of k only. They will be imaginary for real k. Define 
the real valued surface 2 to be the locus of points (k, -j\ (k)) in four- 

A.. 1 

1 

space. Then for a given go = co the intersection of the two surfaces 

o K 


e 


H- o 


/ 5 , 


ml 


rr2 



(6.51) 


y.Po, 



XmPi 


N' ON 


and 
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-jX_ - Gl> o = 0 in four-space is a curve. But this curve in four-space is 

equivalent to a surface in three-space ( k space) which is commonly called a 

sheet of the dispersion surface at co = go . The totality of all these surfaces 

o 

in three-space for all i = 1,2,.. 6 + 4N is called the dispersion surface at 

angular frequency go = co^. However } at a given frequency the plane - go q = 0 

may not intersect ; therefore, there may not be 6 + 4N real sheets to the 

r 

dispersion surface at a given frequency. As go takes all real values, the curve 

in four-space which is the intersection of 2 and the plane -j\. - go = 0 

K . i 

th 1 

generates the i sheet of the dispersion surface in four-space. This is 
equivalent to relabeling the -j\ axis as go. 

Recall that for source free solutions it is required that 


det|0~* 

But det U > 0; hence, all of the sheets of the dispersion surface are asso- 
ciated with the eigenvalues of the characteristic Equation (6.52). Therefore, 
there is a one-to-one correspondence between all sheets of the dispersion sur- 
face S. = \. - Jgo = 0 and the eigenvalues of the characteristic equation. This 
is in direct contrast to the six-vector method in which the factor det U(1 c,go) 
may be zero and contribute to sheets of the dispersion surface. Further 
remember that it was possible for f|~ U(k,Go) = 0, for some i of the six- 
vector method to coincide with a sheet (s) of the dispersion surface. But 
for the 6 + 4N-vector method, U F. = 0 if and only if F. = O since U is 
positive definite. This one-to-one correspondence between the sheets of the 
dispersion surface and the eigenvalues is one of the advantages of the 6 + 4N- 


cjyJ = IdetipTr U,--jcj) = 0 


(6.53) 


vector method. 
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Another property of the dispersion surface can be determined from the real 
valued property of the fields, i.e., both ^ and are real valued vectors. 
The real valued property of the field implies that 


Olk) = 0* l- k) 


(6. 54) 


Therefore , 


det 0(k) — X Ik) 


uj = det jo* 


(-k) — XI k )U =0 (6.55) 


Which implies that 


det [oho - (-X(k)) u] = 0 


(6. 56) 


and finally 


det [ Olk) - (-Xl-k))u] = 0 


(6. 57) 


Hence if X^Ck) is an eigenvalue, there exists a j such that \ (k) = - X^Ok) 

is also an eigenvalue. From this we have that the i^* 1 and j** 1 sheets of the 

dispersion surface are S. = X . (k) - jco = 0 and S . = - X.(-k) - ju> = 0, res- 

11 J i 

pectively. Thus, if either point (k,o>) or (-k,-oo) is on the dispersion surface 
then both points are. 

The 6 +• 4N eigenvectors which are a function of k only also display an 
orthogonality property. This property is 


F: + IJF. = O for \ . 


i j = l,2,"‘ 6 + 4N 


(6. 58) 
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But even for a degenerate system one can define a set of orthogonal eigen- 
vectors. Thus, it will always be assumed that the eigenvectors form a mutually 
orthogonal set, i.e. , 


FfuFj = 0 


i,j= 1,2/ 


• 6 + 4N 


( 6 . 59 ) 


Another property may be surmised by examining the original real system of 
equations. In a lossless system the property of time reversal essentially 
states that if the time axis were suddenly reversed then the system would re- 
trace its path of operation. That is, time reversal is expected if the system 
is invariant under the following transformation. 


t — 

► -t 

t — 

— g 

— 

ft 

— 


— 

— - u„ 

3o — 

-B 0 


That time reversal is expected for our plasma model is easily verified. Thus, 
it is anticipated that time reversal should manifest itself in some manner in 
the eigenvalues, eigenvectors and the dispersion surface. If the Fourier 
transform of a function f(r,t,B Q ) is g<k,co,B o ) , then the Fourier transform of 
the time reversed function + f(r,-t,-B Q ) is + g(k,-a),-B Q ) . Therefore, for a 

reversed characteristic field is as follows. 


characteristic field the time 
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TABLE II 


E j ( k , X j , B 0 ) 


Hi (k.X,. B 0 ) 


v„i( k , X , , B 0 ) 

u n j ( k , X j } Bq ) 


Characteristic Field 


Ei (k, Xj ,-B 0 ) 

-H f (k, X j -B 0 ) 

-x, 

-V nf (k,X,-,-B 0 ) 

“ni^, Xj,-B 0 ) 

Time Reversed Characteristic Field 


It is easily verified that if the characteristic field satisfies the oigenv.i I o, 
equation, the lime reversed characteristic field does also. Hence, if the 
characteristic fields are known for all i such that -j\. > 0 , then the oilier 
half is known by the property of lime reversal. With respect to the dispersioi 
surface, time reversal implies that if point (k,u>) is on the dispersion surUe. 
then (k,-») is aJ so. Time reversal together with the real valued proper! y 
"uplus that if point (k,«) is on the dispersion surface tlien point (-k,o>) is 
also. Therefore, the three -space dispersion surface at a given o; has the 
symmetry of reflection through the origin. In addition, time reversal and the 
real valued property implies that if any one point of the set (k ,o>) , (lT,-io), 
(-k, w ) and (-k,-*o) is on the dispersion surface then all points are. 

Since the group velocity of the i th characteristic wave is 

V;Si Vj-Ui-io)) • 


v gi 5 


did 




= -iVtx, 


( 6 . 61 ) 
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i lu-n the group velocity of the time reversed characteristic wave is 

Vfi-x, -jw) 


d 


(-X j -j OJ) 


(6.62) 


dcu 

= iV T x, 

* ~ v gi 

Then lore, if a characteristic wave at (k,co) satisfies the radiation condition, 
then the time reversed characteristic wave at (k,-c o) does not. And 
will satisfy the radiation condition but (-k,co) will not. Thus, when the exact 
tc»t al field is sought, the sum (integral) of characteristic waves must be 
taken in such a manner that points (k,co) and (-k,“u>), which satisfy the radia 
tion condition, are included and points (k,-co) and (-k,co) are excluded. This 
may mean that the contour of integration must be chosen appropriately. 

6 6 Case of a Compressible Plasma with N = 1 Species of Charged Partic I o s 
(10-Vector Method) 

For the case of a mobile electron with an isotropic pressure, the 10 th 
order determinant al equation reduces to 


-\€ 0 -jkx 


jkx -X Mo 


* 

0 


-P r 


q 0 - ^qB 0 X 

rT 




-xyp 0 


(6.63) 


= 0 


After considerable algebraic manipulations it is found that the equation becomes 
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(yP 0 )(^m 


x'°+ |( 2+ S) y 2 + 3 cUn + cJh} ^ + {u + 28)y 4 + (4 + 8 ) CUNY 


+2cuHy 2 + S(aJh£ 2+ w h w n + 3 Ci)n| ^ + |sy + 28cJny C 
+ j^(l + 28)0U|y|+ cu^ y + (2 + 8)cUjg + GJnCiJivijy + gUh^nC + 0Dn^^ 

{w 2 ^C 2 y 2 + Sw 2 £ 2 y 4 }x 2 J =0 


where OJfj - /^m , 8 - 7 P 0 /^°m c 

and y = Ck=(£,77 1 £). 

The eigenvectors corresponding to the nonzero eigenvalues are 

E = P m (\P q )~'(\ Z + X cD H x + 8 y y T ) V 
H = jPq C ( XO>^) l_ y X ( X + CU H X)V 
V = V 

u = j (Xc) '"y T V 


(6.64) 


(6.65) 


( 6 . 66 ) 


(6.67) 


( 6 . 68 ) 
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or explicitly 

1 


E x = (A> m /x^ q )[(X 2 + se 2 )v x+ (-Xcu H + %£rj)\f y + 8£C v z] 

(6.69) 


Ey = (^ m /xP q )[(XW H +S£t 7 )v x + (X 2 *87 7 2 )V y + S7 ? tV z ] 

(6.70) 


Ez » (/ > m/^q)[se^V x + Sl7^V y MX 2 + 8C 2 )V 2 ] 

(6.71) 


h x = (j/> q c/xo)N)[-C^HV x -£xv y + r/xv z ] 

(6.72) 


Hy - ( j Pq C /x OJjj ) jjj XV x - £ V y - £ XV Z ] 

(6.73) 


H z = U^qc/x(j^)[((aj H -T]X)v x * (7]a; H +CX)v y l 

;i 

(6.74) 


V x = X 2 {x 5 - |\i - SK 17 2 +£ 2 )- 2 ( y 2 + Od^) j X 3 + OJ H t, Tf X 2 

-(y 2 -aJN}|(i-a)(97 2 4-^)-(y 2 + ^)jx + cu H CT?(y 2+ 4 js| 

(6.75) 


v y = x 2 {-aj H x 4 >li-a)T7CX 3 -cu H [2y 2 -77 2 -(i-8)C 2 + ajNjX 2 

■*-(l-8''7^y 2 +cDN)>“W H [(y 2 - < »7 2 )(y 2 -*-^S-0-8)y 2 C 2 ]} 

(3.76) 


v z = £x{(i~ 8 )£X 4 + 8 o)h t }\ 3 + eju-SHy 2 + u .? 2 j " 

r "1 r *> 

+ 1 ; H T) jj$y 2 «• tofi \ X + Wh Tv 2 ^ 

(6.77) 


u = jlXcf 1 £V X + v Vy + £V 7 j 

(6.78) 
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The remaining two eigenvectors corresponding to the eigenvalues which are 
identically zero are 


Fp = 0‘, k, 0,0 


(6.79) 


k,0,0, )P q 


A P o_ 


(6.80) 


The actual process of computing the patterns at a given frequency go and 
distance r is not difficult. Only simple algebraic and differentiation opera- 
tions are involved. However, even for this simplest warm plasma model either 
a computer or fortitude must be used because of the complexity of the alge- 
braic quantities. The foJ lowing is a brief outline of the steps involved :i n 
computing the "pattern." 

(1) Determination of the dispersion surface. From the determinant al 
Equation (6.64) it is seen that the dispersion surface in k-space is a surface 
of revolution and has at most three real sheets. Therefore, the dispersion 
surface may be found by obtaining the six roots of y with 0 = cos ^ (t>/y) as 

a parameter ranging from zero to 277 . 

(2) Determination of stationary points for a given r. The stationary 
points are those points in which the plane tangent to the dispersion surface 


is perpendicular to r. This is equivalent to points on the dispersion surface 

A d f 1 

such that Vj^(k*r) = 0, which reduces to — [k(0) cos (0 - a)J = 0, (a is 

A A 

the angle between z and r) since the surface is a surface of revolution. Hence, 


A A 

the^ stationary points lie in the plane determined by r and z. The stationary 
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points that satisfy the radiation condition must also satisfy the condition 

A 

k* r > 0. The stationary points may be found either by use of a computer or 
graphically. 

(3) Determination of the Gaussian radius of curvature. The Gaussian 
curvature is defined as the product of the two principal curvatures 
Pi P 2 • Since 2 is a surface of revolution p^ (meridian curvature) and p^ 
(parallel curvature) are given by 

P\ = (k 2 ♦ 2 k /2 - k k") ( k 2 ♦ k' 2 f 3/2 (6. si) 


-i 


Po=k sin a esc 9 


(6.82) 


where the prime denotes the derivative with respect to 0. Point k on J is a 
stationary point; hence. 


k ‘ cos (a-0) = - (k') 1 Sin (a- 0) = (k 2 + k /2 ) 1/2 


(6.83) 


From Equation (6.83), one finds 


= (k 2 + 2k /2 -kk ,, )(k 2 + k ,2 f l 

and 


^1 = k COS (Of- 6) ^ 
Then it immediately follows that 


(6.84) 


(6.85) 


ftsk-'kl'if 
0 d 9 


( 6 . 86 ) 
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where 


sin R c sc Cl sec (Q i- 8) 


(6.87) 


(4) Normalization of the eigenvectors. The normalization process 


Re 


EXH* + yp 0 u V* 


= I 


( 6 . 88 ) 


as a general function is a rather formidable task, however, is simple pointwise 
on a computer. 

Needless to say, tbe remainder of the pattern calculations are rather 
trivial . 

6.7 Case of a Cold Plasma with N - l 

Although the cold plasma is a particular case of the previous warm plasma 

and may be found accordingly, drastic simplifications make it quite feasible 

to find the radiation "patterns" without the use of a high-speed computer. 

2 

First the determinant al equation becomes quadratic in k and the usual Appleton- 

Hartree equation results. The dispersion surface then can be found by a 

22 

graphical method of Deschamps and Weeks, Also both the stationary points 
and the radii of curvature can be found graphically. Furthermore, the polari- 
zation ellipse of Fb has one axis m the plane (k,lO and the other perpendicu- 
23 24 

lar to it. ’ If the axial ratio is tan (3, the angle (3 is simply related to 

2 2 

0 and to the usual ionospheric parameters X = oi./gj and Y = co__/gj by 

N n 

ten 2/3 = 2 ccs 9 csc*"0(l~X)Y 


(6.89) 
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The two solutions give the polarization ratios to the "ordinary" and "extra- 
ordinary" waves. For the usual definition of the "ordinary" wave (the sheet 
of the dispersion surface for which k = k ~\/l-X when "k-H =0) the major 
axis of the H-ellipse is perpendicular to the plane (k,lI o ). The major axis 
of the H-ellipse for the "extraordinary" wave lies in the plane (ik H ). The 
normalization requires that 


|Hj| 2 = (k /0JfJ. Q ) COS lOf-0) (6.90) 

The vector E^ follows form Maxwell's equations 

E| = - (I/Cl)) 6 'kxH. ( 6 . 91 ) 

The above brief description along with Equation (6.31) should prove sufficient 
to find the radiation field of an arbitrary antenna in a cold plasma. 

The principal advantages of this method of obtaining the radiation fields 
of an arbitrary antenna, over that of previous authors are: 

(1) A high-speed computer is not necessary. 

(2) Functions of an indeterminant form at the principal directions such 
as in Bunkin' s and Keuhl's solutions are not involved. 

(3) All terms of the solution and f e normalization are easy to 
interpret physically. 

It would be interesting to display the dispersive property of the medium by 
showing a diagram of the four-space (k,-jX.) dispersion surface, however, such 
is not possible. Fortunately since the k-space dispersion surface is a surface 
of revolution, the dispersive property of the medium may be illustrated by a 
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diagram of (k ,k , -jX). Figures (6.2), (6.3), (6.4) and (6.5) show the 
p z 

sheets of one quadrant of such a surface in which co H ./w N = 2. The sheets of 
the dispersion surface in the other quadrants are mirror images in the quad- 
rant planes. Since there are nine eigenvalues, the remaining sheet is -jX = 0. 
With the dispersion surface represented in this manner, there are only two 
topologically different dispersion surfaces, I 1 = Notice also that the 

sheets touch only at -jX = 0, + u> . Because of the complexity of sheet 2) 

N A. g 

in Figure (6.3), perhaps a verbal description is in order. There are four 
values of -jX that correspond to important points on the sheet. The sheet 


intersects the -jX axis at 


1 2 . 2A 

-JX = 2 [C^h + 4 V " “h 


which is called cutoff. 


At -jX = w , 2. intersects the light cone with the portion of the surface for 
N X 2 

-JX < "inside” and for -jX > oo^ "outside 9 ' the light cone. -jX - co^ is an 

asymptote for the intersection of S with the (k , -jX) plane. And 

2 Z 

O 9 r 

-jX - (oo^ + c^ N ) 2 is an asymptote for the intersection of 2^ with the (k^, ~ X) 
plane. In the neighborhood of -jX = 2^ changes most rapidly. There is 

a "plateau" at -JX = oo^ which is widest in the direction and disappears in 


the k direction. 


P 

The intersection of planes -jX = go with the dispersion surface is what is 

called a CMA diagram (Clemmow, Mulally, Allis). They are represented by the 

dotted lines in the figures. The intersection of a plane passing through the 

25 

-jX axis with the dispersion surface yields a Stringer diagram, which is 
represented by an alternate dot-dash line. The curves corresponding to either 
side of the removed sector is a Stringer diagram. 

One of the advantages of the dispersion surface represented in this manner 
is that it is emphasized that there are more than one "ordinary” and extra- 
ordinary" sheets. Sheets 2 an d 2 correspond to the "ordinary' sheets 

X i x 3 







141 


while £ and £ correspond to the "extraordinary" sheets. Also sheets 
X 2 X 4 

£ and £ are of left hand polarization while £. is of right hand polari- 
X 1 X 4 o 

zation. Furthermore, the portion of £ "inside" the light cone (-j\ < col.) 

2 „ „ 

corresponds to right hand polarization while the portion "outside" (-j\ > V 

corresponds to left hand polarization. 

Briefly, let us consider the major differences in the dispersion surfaces 

of the warm and cold plasmas with one mobile charged particle species. First 

the warm plasma system matrix can be expressed as a bordered cold plasma system 

matrix. Hence, the ten warm plasma eigenvalues will alternate with the nine 

cold plasma eigenvalues. Likewise the sheets of the warm plasma dispersion 

surface will alternate with the sheets of the cold plasma dispersion surface. 

The sheets of the warm plasma dispersion surface do not touch at -j\ = u> N but 

may touch elsewhere. And the warm plasma sheet corresponding to £ is not 

1 f 2 2 i 2 2 i 

limited to values of ~j\ between — (u> H + 40^) 5 and + but 

extends to infinity in the -j\ direction. That is, this sheet for largo 


k, -jX becomes large. 
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7, CONCLUSIONS 

A phenomenological approach has been used to investigate some properties 
of linear passive media. Properties of fields that can propagate in linear 
passive media were postulated and from this properties of the media through 
the Four ; or transform of the constitutive matrix were deduced. The concept of 
o positive real condition on the constitutive relationship for linear passive 
media was introduced and some of its implications were considered. Also, the 
concept of causality,! which is more fundamental than the group velocity concept 
and which is necessary for realizable media^was considered, pari.ici ,r r!v r for 
t he e as e of i s o t r op :i c me d > n . 

A general formulation of the spectrum of characteristic waver: in loss 
linear passive media has been made. The fields due to an arbitrary source can 
be separated into components parallel to the characteristic waves by using an 
orthogonality condition for the characteristic waves of the medium. The coin 
porents of the source field are dependent only upon the portion of the source 
parallel to their characteristic field and to their own sheet (s) of the dis 
pci's ion surface. The theory has been applied to the problems of longitudinal 
and transverse electric dipoles in a general time-dispersive uniaxial medium 
and to an electric dipole in an isotropic compressible plasma with N species 
of charged particles to obtain exact solutions. 

Finally^ the radiation field of an arbitrary source in a lossless linear 
passive medium has been obtained by using the spectral decomposition of the 
fields and the stationary phase method. It is shown that by normalizing the 
length of the Total Poynting vector (electromagnetic plus medium) to unity 
for each characteristic field, a concise and physically interpretable 
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expression for the source fields is obtained. The radiation field for a time 
harmonic source was found to depend upon the Gaussian radius of curvature, the 
reaction of the source with the normalized characteristic field, and the char- 
acteristic field at each stationary point on the dispersion surface. These 
results have then been applied to an anisotropic compressible plasma and to a 


magneto-ionic plasma. 
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